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THE SUN SOURCE COSMIC RAYS 
INTERMEDIATE 


KatzMan 


ABSTRACT 


The cosmic ray intensity measured with extremely narrow-angle tele- 
scope, steradians, and with inches lead absorber for the period 
January 1955 December 1958 shows increase 20%. This increase 
attributed particles coming from the sun. shown that the change 
hour maximum the first and second harmonics can explained change 
the relative importance the impact zones. This phenomenon attributed 
change the number and polarity sunspots. 


INTRODUCTION 


The long-term change the intensity cosmic rays has been shown 
Forbush (1954), Fenton a/. (1958), and others relate inversely with solar 
activity. The activity has pronounced influence the intensity particles 
| 


low energies monitored neutron monitors, ionization chambers, and 
Geiger counter telescopes. The variation cosmic rays high energies 
studied means extensive air showers but such experiments have not 
operated continuously for complete solar cycle. avoid confusion will 
call those energies particles between that monitored 
Geiger counter telescopes with the conventional lead absorber 
and those energies particles studied air showers. The experiment 
Rau (1939), who lowered ionization chambers into Lake Constance depth 
meters during the solar maximum will fall into this category. 
The very interesting fact reported, though unexplained, was that the second 
harmonic the cosmic ray intensity showed positive correlation with the 
second harmonic atmospheric pressure. 

The activity the sun known pass from minimum maximum 
period approximately years and measured terms the number 
sunspots, the magnetic polarity which has period approximately 
(The Sun 1953). Thambyahpillai and Elliot (1953) have suggested that 
the change hour maximum the daily cosmic ray intensity may 
related the 22-year sunspot cycle. Since the sun active body and 
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continuously emitting energy one form another, reasonable 
consider that the sun under certain conditions does emit particles. When the 
activity low ebb, measured solar cycle, the particles will 
low energy and will either not leave the influence the sun those that 
reach the outer atmosphere the earth will not have sufficient energy 
that their progeny will reach the surface the earth. this study 
consider the solar minimum activity that period when particles from the 
sun not have sufficient energy influence monitors the earth and the 
intensity cosmic rays during this period that particles that come from 
outside the solar system. 

The magnetic field that surrounds the earth more that dipole 
some distance from the earth than comparatively small distances from 
the earth. The field that cosmic ray particle encounters not that the 
earth’s field only (Rose al. 1956). has been shown Dwight (1950), 
Firor (1954), Liist (1957, 1958) that the magnetic field considered 
that dipole then certain hours the day (impact zones) will show 
greater intensity cosmic rays than the other hours particles are arriving 
from the vicinity the sun. They found 0300, 0900, and 2000 the import- 
ant hours. radical change the magnetic field that the cosmic ray particles 
see they approach the earth will change the relative importance the 
impact zones. 

explain the sudden decreases intensity cosmic rays known 
Forbush decreases, which last for number days, and the long-term 
cycle considered created the sun’s activity that 
bars, either electrical magnetic means, particles low energy from 
reaching the earth (Simpson 1957, 1958; Morrison 1956; Alfven 1954; and 
Parker 1956). This barrier must maintained and intensified with solar 
activity explain the cycle cosmic ray intensity. explain the 
Forbush decrease the barrier must increased suddenly for short period. 
field force can exist unless associated with matter one form 
another. accept the hypothesis that the sun emits particles all times 
and that barrier created this activity then the following can 
expected: 

(1) decrease low-energy particles monitored neutron monitors, 
the barrier stops more low-energy particles from reaching the earth than 
the sun emits; 

(2) net decrease particles monitored Geiger counter telescopes 
ionization chambers with lead absorber the barrier intensity 
increases; 

(3) increase intensity short duration before Forbush decrease, 
the intensity the barrier being increased after the particles pass the earth; 

(4) possible increase intensity monitored telescopes not sensitive 
particles energies below the intermediate range. 


EXPERIMENTAL 
narrow-angle telescope steradians) with inches lead 
absorber has been operating continuously since December 1954 (Katzman 
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1958). Triple coincidences are recorded 20-pen Esterline Angus tape 
recorder and scalers. The data, for the period January 1956 Decem- 
ber 1958, were analyzed and are presented the form curves (Figs. 1-3). 
The data were not corrected for either temperature pressure and are shown 
12-month running averages. 


RESULTS 

The 2800 Mc/sec radio wave, emitted from the sun, used measure 
solar activity (Covington 1957). The curves Fig. are 12-month running 
averages percentage the 1955 average value. Curve plot the 
radiation intensity frequency 2800 Mc/sec,? curve the critical 
frequency the and curve the intensity cosmic rays 
monitored the narrow-angle telescope. 12-month period data 
referred the first day the seventh month that period, thus July 


flux 2800 Mc/sec sec (bandwidth Data supplied 
Mr. Covington, National Research Council, Ottawa, Ontario, Canada. 

supplied Defence Research Telecommunications Establishment, Defence Research 
Board, Ottawa, Ontario, Canada. 
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1955 refers the 12-month period January December 1955. There 
1956, and curve September 1956. clear that about Sep- 
tember 1956 radical change took place the conditions existing the 
sun-earth system. The rate change the solar flux then decreased. The 
critical frequency the layer practically stopped increasing and the 
cosmic ray intensity remained constant until March 1957, then dropped 
sharply August 1957, and later rose slightly. 

The cosmic ray intensity data are plotted hourly basis Figs. 
and Fig. each hour the day plotted separate curve, the hour 
for 12-month running averages taken percentage the mean for that 
12-month period. The 9-10 hour only 80% the annual mean for 1955 
and rises 105% for 1958, whereas the 15-16 hour 116% the 1955 
mean and drops 86% the 1958 mean. The daily intensity plctted 
percentage the annual mean Fig. again each curve represents year 
data centered about the day indicated. The intensity about the 0300, 
the 0900, and the 1600 and 2400 hours above the annual mean. The 
pattern curve No. maintained curve No. when there marked 
change and the pattern changes again curve No. 25. 

can readily seen from the curves Fig. that the hours maximum 
for the first and second harmonics obtained from curve No. would not 
the same that curve No. 37. Since the 15-16 hour intensity con- 
sistently high the 20th period (referring now Fig. and then reduces 
uniformly low value, the month when the value reached 100% the 
annual mean was chosen the dividing month. The data were thus divided 
into the two periods: January 1955 August 1957 and September 1957 
December 1958. harmonic analysis for the first and second components 
the diurnal variation for these two periods shown Fig. The hour 
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dials the diurnal and semidiurnal components the atmospheric 
pressure, Ottawa, Canada, and the cosmic ray intensity for the two periods, and 
January 1955 August 1957 and September 1957 December 1958. 


maximum the first harmonic hours minutes for the first period 
and hours minutes for the second period The second harmonic 
also shows large change, the times maximum being hours minutes 
and hours minutes for and was mentioned above that these data 
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were not corrected for barometric pressure. harmonic analysis the diurnal 
variation barometric pressure does not show any appreciable change over 
the period. Temperature effects the atmosphere though they have seasonal 
effect could not account for this long period change. clear, therefore, that 
this change phase not atmospheric origin. fact correlation with 
atmospheric pressure would show change from negative positive. The 
pressure vectors are also shown Fig. 


DISCUSSION 


The increase ionization the layer cannot attributed only 
increase electromagnetic energy emitted the sun, since the percentage 
the increase over the earth not uniform (Katzman 1958). Therefore, 
some the increase ionization must due increase the number 
particles arriving the earth. Since has been shown others that the 
sun’s activity bars low-energy particles from arriving the earth then this 
increase particles must coming from within the system. 

The intensity cosmic rays intermediate energy remained constant 
February 1956 and increased 40% September 1956, remained 
steady till April 1957 when dropped level 20% above the 1955 annual 
mean. This increase cannot attributed increase the cosmic ray 
intensity outside the barrier created the sun. The increase must, therefore, 
attributed some mechanism within the barrier. Any mechanism which 
would increase the intensity intermediate energy particles would also 
increase the intensity the low-energy particles. This can also seen 
the nucleon component Ottawa from June December 1956 Fig. 
Fenton al. (1958). The sunspot number continued increase during 
this period whereas the nucleon component showed increase instead 
the expected decrease. Hence, since neither the sunspot number nor the solar 
flux showed decrease the increase the cosmic ray intensity inter- 
mediate energies for this period must attributed some mechanism 
the sun. 

Under certain conditions possible for the sun’s activity the cause 
increase cosmic ray particles energies above Bev was evidenced 
the sudden increase February 1956 Huancayo (Simpson 1958). 
study this behavior the intensity the nucleon and meson components 
Ottawa shown percentage the average intensity the first days 
January 1956, for the months January, February, and March (Fig. 5). 
There reason believe that the weather conditions for similar months 
other years will drastically different that the intensity would 
radically different because the weather; thus similar periods for the years 
1957 and 1958 are also shown for comparison. These curves are very inter- 
esting. Before the ‘‘flare the meson component remained above 
100% for about days while the nucleon component after small decrease 
slowly increased intensity peaked value equal the meson com- 
ponent just before the Forbush-type decrease. The decrease the nucleon 
component preceded the decrease the meson component and the reverse 
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Fic. The intensity Ottawa the nucleon and meson components, corrected for 
pressure, percentage the first days January 1956. 


obvious for the return intensity. The following decrease 
March not marked for the meson component but the increase before the 
decrease clear the meson component. The intense decrease January 
1957 preceded increase both the nucleon and meson components 
centered about January. the intensity the barrier increasing 
that particles higher energy are being stopped from entering the 
system, then with time decrease particles energies between the low 
energy the nucleon component and those intermediate energies main- 
tained. This can seen for the years 1957 and 1958, Figs. and when 
the intensity the meson component dropped relatively more than the 
nucleon component. Since the intensity the cosmic rays intermediate 
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Fic. The intensity Ottawa the nucleon and meson components, corrected for 
pressure, percentage the first days January 1956 for the months October, November, 
and December 1958. 


energies increased about 20% while the intensity particles lower 
energies decreased proposed that the sun some means yet un- 
explainable causes barrier, suggested others, created that stops 
cosmic ray particles lower energies from outside the system from reaching 
the earth, while the same time emitting particles energies inter- 
mediate energies. Thus monitors that are sensitive lower energies show 
decrease whereas monitors that require higher energies show increase. 
The experiment Rau (1939) Lake Constance was performed 
approximately the same geographic latitude Ottawa. The hour maximum 
for the second harmonic atmospheric pressure would, therefore, about 
the same local time for the two stations. Since obtained positive corre- 
lation with pressure for cosmic rays intermediate energies (the hour 
maximum for the second harmonic was 1000 hours) similarly reported 
this paper reasonable conclude that the cause that increased the 
intensity the morning hours and decreased the intensity the afternoon 
hours was the same for both the solar maximum and 1957-58. 
The condition the sun that known the same for these two maxima 
the magnetic polarity the sunspots. This indicates that the magnetic 
intensity and polarity the sunspots move from the sun towards the earth 
carried particles emitted from the sun. This influence would greatest 
the sun’s maximum activity since then that the greatest number 
spots exist the sun. the sun’s activity decreases the number spots 
decrease and the influence the polarity the spots decreases till eventually 
the earth’s field becomes dominant. Rerurkar (1957) with the aid 
(1954) theory frozen-in magnetic field ionized beam emitted 
the sun has shown that the hour maximum will shift from the afternoon 
morning hours when the direction the magnetic field rotates through 
180°. This clearly shown Figs. and where the predominance the 
16th hour intensity shifts with increase solar activity the 11th hour. 
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CONCLUSION 


From the data presented concluded that the sun emits cosmic ray 
particles whose energy depends the intensity the solar activity. The 
magnetic field the sunspots carried towards the earth the emitted 
particles. The change number and polarity the spots will influence the 
impact zones the surface the earth cosmic ray particles and will cause 
apparent change correlation with atmospheric pressure from negative 
positive. This will also cause change apparent anisotropy with change 
solar activity. 
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RESOLVING POWER THREE ANTENNA PATTERNS 
DERIVED FROM THE SAME APERTURE! 


ABSTRACT 


Three antenna patterns are derived from the same linear aperture and may 
described terms angular spectrum spatial frequencies ranging from 
zero common frequency. The band passes according the shape 
the spectrum are rectangular, triangular, and cosinusoidal for the three patterns, 
and give resolving powers respectively 1.33, 1.00, and 1.05, terms the 
period. The rectangular band pass gives rise the optimum antenna pattern and 
allows the Fourier components source from zero cutoff frequency 
received with equal intensity and zero phase shift. Scanning curves two equally 
intense point sources and uniformly bright line are investigated. 


INTRODUCTION 


radio astronomy, the abundance and closeness radio-emissive objects 
the sky has necessitated the construction various types antennas with 
large dimensions, effort isolate and study the size individual 
objects. Although the apertures these antennas, now use under con- 
struction, are measured terms several hundreds meters, minute 
arc about the narrowest beam width which can realized. There are objects 
the sky which are much wider than this, but other sources such certain 
sunspots and radio stars may have widths the same order. Thus, the latter 
case, necessary attempt obtain true picture the structure the 
object under examination, least some estimate its structure, the face 
serious limitation antenna resolving power. This problem has two 
aspects: first, that utilizing given physical area obtain antenna 
pattern with the sharpest main lobe; and second, that interpreting the 
departures the observed scanning curves from the known antenna pattern 
terms object structure. Both these aspects will discussed, since they 
are interrelated. The simplest formulation the problem, exemplified the 
fan-shaped antenna patterns, will developed and will have the advantage 
introducing the simplicity single variable. The resolving capabilities 
three different antenna patterns, each derived from the same over-all 
aperture, will found. Two these patterns have been realized already with 
slotted waveguide array (Covington and Broten 1954; Covington and Broten 
1957), while the third remains developed. 


THE OBSERVED BRIGHTNESS PROFILE OBJECT 


The general equation relating the observed temperature distribution pro- 
duced scanning distant object with fan-shaped radiation pattern 
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the true temperature distribution this object may written for the single 
dimension under consideration 


where the equivalent temperature the antenna radiation resistance 
when the maximum the antenna pattern angular distance from some 
reference point, 

the normalized power pattern the antenna function 
the angular distance from the maximum the antenna pattern, 

S(x) the temperature brightness the source, function the 
angle from the reference point. 

This equation shows that the observed brightness profile obtained from 
scanning the object with the antenna pattern differs from the true one 
result averaging process introduced the antenna pattern itself. 
alternative relation between the Fourier transforms the spatial frequency 


variable given the equation 

(2) R(s) S(s), 

where R(s), A(s), and S(s) are the Fourier transforms for the three functions 
defined for eq. (1); for example, 


The variable the number angular wave crests described per unit angle 

The use the concept spatial frequency suggests the analogy that the 
antenna may regarded filter and that the spectrum the antenna 
pattern may spoken the antenna band pass. For antenna whose 
aperture large, measured number wavelengths, the band pass has two 
features: 

(i) zero for all values greater than cutoff frequency given the 
relationship 


(4) 


where the length the basic aperture under consideration, and the 
wavelength observation; 

(ii) takes definite shape the region below cutoff, attenuating some 
frequencies relation others. 
the process making observation, these two characteristics the 
antenna band pass make their imprint the spectrum the object and 
account for all the distortion introduced into the true profile the object. 
terms frequencies given eq. (2), this amounts the complete removal 
all frequencies which are greater than cutoff frequency and, general, 
changes the amplitude and phase all the lower frequencies which are 
accepted the antenna. 
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The solution eqs. (1) (2) for the source function with known antenna 
pattern which has finite cutoff frequency has been called the principal solution, 
since its spectrum the same the spectrum the true distribution for all 
the frequencies below cutoff (Bracewell and Roberts 1954; Arsac 1955). This 
function may not good approximation the source function, since the 
higher frequencies which were completely suppressed the antenna are 
missing it. However, the Fourier approximation with finite number 
terms and the best that can obtained unless there some other 
edge the source available. obtain the principal solution from observed 
curve taken with antenna which attenuates some frequencies more than 
others, process successive approximation must undertaken that will 
restore all the frequencies below cutoff their original values. This restoration 
process will limited the presence errors, such receiver noise, since 
the restoration increases the amplitude given frequency regardless its 
origin. The principal solution may also obtained directly when the antenna 
has constant band pass from zero the cutoff frequency. this case, the 
power pattern represented the function (sin and 
may approximated very closely practice (Covington and Broten 1957). 
restoration need done with the observed results and this sense 
that this particular antenna pattern with its relatively large side lobes 
optimum. 


CLASSIFICATION ANTENNA PATTERNS MEANS BAND PASS 


The spatial frequency spectrum the antenna pattern provides alterna- 
tive means describing the directional properties antenna and will 
used describe some the various types antennas which have been 
developed within the past few years for radio astronomy. For the purpose 
evaluation, the antennas considered are all characterized the collection 
energy from the same aperture, but are physically different the method 
using this aperture. Fig. some these configurations are shown together 


ANTENNA TYPE CONFIGURATION BAND PASS 

AREA UNIFORMLY ILLUMINATED (a) | 

TOTAL POWER DETECTION | as 
Se 
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° Se 
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Fic. Various types antennas, their physical configurations, and the corresponding 
angular spatial spectra. 
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with the corresponding band passes. all cases, with the exception the 
cross, the band passes have common cutoff frequency. the construction 
the successively larger and larger antennas which are required some 
radio astronomy studies, will readily appreciated that, some physical 
size, further fabrication will possible only through the use individual 
elements, each erected firm foundation and connected together some 
manner such shown the configurations Fig. The important question 
power gain, the methods making the primary elements, and the directing 
the beam will not considered here. 

The basic area available for the collection energy from the incident wave 
shown rectangle Fig. and with suitable antenna design may 
made into single collecting area. When the voltages contributed all points 
the aperture are equal intensity, the power delivered the receiver will 
depend upon the direction arrival the wave given the receiving 
power pattern defined equation 


(5) 


where the variables have been defined previously. The corresponding band 
pass the spatial frequency spectrum, Fig. triangular function with 
maximum amplitude zero frequency and minimum the cutoff frequency 
This cutoff frequency the maximum which may obtained from the 
given length the basic area, when greater than unity, and related 
eq. (4). All the other antenna patterns which may obtained 
placing different emphasis the manner collecting voltages from various 
points along the aperture have band-pass spectra which are contained within 
this frequency. 

The collecting areas shown Figs. (d), (e), and belong various 
types phase-sensitive interferometers, all having the same antenna pattern 
the short dimension but with various types patterns the long dimension. 
Again uniform collection voltages from the incident wave assumed for 
each the antenna elements. the phase-sensitive interferometer, the 
resultant power pattern the product the voltage patterns the two 
interferometer elements multiplied the usual interference term which 
related the distance between the phase centers and obtained through 
electrical process. This based upon slow, continuous shifting the phase 
angle between the two radio-frequency voltages picked the interfero- 
meter elements that these voltages are successively phase and then out 
phase. This leads the formation two antenna patterns the receiver 
output, one associated with the mean power and the other with the fluctuating 
power. The desired antenna pattern associated with the fluctuating power 
isolated means low-frequency phase-sensitive detector operating 
synchronously with the slowly rotating radio-frequency phase shifter. this 
preliminary discussion, the error arising from the fact that necessary 
break shorten one the slit apertures the interferometer elements will 
not considered, since its effect can made small and compensated for. 

The simple interferometer with the maximum separation permitted the 
basic area under consideration, Fig. has band pass consisting single 


2 
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line the cutoff frequency where the band pass the uniformly illuminated 
aperture, Fig. goes zero. The difference response this frequency 
gives rise the statement that the interferometer pattern has twice the 
resolution the uniformly illuminated array for, although multilobed, the 
lobes are much narrower. This pattern is, course, only useful when 
restricted the primary patterns the interferometer elements region 
known contain single source. 

one form the variable spacing interferometer, Fig. 1(c), one the 
antennas mounted track that many stations the length the 
rectangle may occupied turn that interferometer element. These 
spacings produce antenna patterns having spatial frequencies from zero 
the maximum defined the simple fixed interferometer, Fig. this 
manner the various spatial frequencies the source may observed without 
distortion and the resultant curves may combined appropriate calcula- 
tions into single curve representing the profile the object. The usefulness 
this method naturally depends upon the constancy the source within the 
interval required take the many scanning curves. 

Configurations 1(d), (e), and make use the same basic physical space 
the others but produce antenna patterns with single main lobe. Fig. 1(d), 
the collecting area occupied two slit apertures the form which 
are operated phase-sensitive interferometer. The antenna pattern this 
case given (sin The spatial band pass rectangu- 
lar, all frequencies being accepted uniformly from zero cutoff. The 
advantage over the variable spacing interferometer that the information may 
obtained single scan the object relatively short duration, thus 
eliminating the necessity making subsequent calculations. 

When two slit apertures the form cross, Fig. 1(e), occupy the basic 
area, the rectangular feature the band pass maintained but the frequency 
response reduced factor two. This reduction resolution arises since 
the phase centers the two interferometer elements overlap, and the inter- 
ference term vanishes completely. 

The compound interferometer, Fig. shows the utilization the area 
under consideration three elementary collecting areas. The antenna pattern 
this case the same that given the and the band pass uniform 
with the same cutoff frequency. physical advantage this type construc- 
tion that the single long element only one-half the basic length. 
similar manner the principle compound interferometer may extended 
four-element grating antenna combined with continuous slit aperture 
length equal the grating gap order produce single-lobed antenna 
pattern (Barber 


THREE SPECIAL ANTENNA PATTERNS 
the discussion antennas which use the same length given rectangle 
for their construction two types band pass with the same cutoff have 
appeared: one which the higher spatial frequencies are linearly attenuated 
form triangular band pass, and second which all accepted frequencies 
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appear equal strength form rectangular shape. This second band pass, 
for reasons previously stated, has been called the optimum for the length under 
consideration. 

great many other band-pass shapes would appear possible, which 
the same cutoff frequency maintained but with various emphasis placed 
upon the relative strengths the higher frequencies. third band pass which 
follows cosine variation from maximum response zero spatial frequency 
zero response cutoff will considered providing shape intermediate 
between the first two mentioned. This band pass and the corresponding antenna 
pattern may conveniently realized with operations based either upon 
antenna with optimum pattern upon the resultant scan obtained with 
such antenna (see Appendix). The first process involves providing for the 
simultaneous reception with two optimum antenna patterns displaced with 
respect each other that the maximum one falls upon the first zero 
the second (an angular shift equal one-half the cutoff period, 
defined later), while the second may conveniently realized the process 
smoothing the scanning curve obtained with single optimum antenna pattern. 
This involves the shifting the observed scanning curve the angle which 
separates the two antenna patterns and its addition the original curve. The 
resultant curve would equivalent that obtained with the composite 
antenna pattern. 

The three antenna patterns corresponding the three types band pass 
are shown Fig. Antenna patterns (2) and (3) would appear more 
desirable than (1) since they have much smaller side lobes. However, the 
larger side lobes (1) are related the uniform band pass and are needed 
record the unmodified high-frequency detail the object being examined. This 
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Fic. Three antenna patterns derived from the same aperture and their corresponding 
spatial spectra. 
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antenna system with the rectangular band pass enables the Fourier components 
source from zero cutoff frequency received with equal intensity 
and zero phase shift, and produces scanning brightness profiles which may 
regarded the Fourier sum approximation with finite number spatial 
frequencies. The use this pattern thus eliminates the need for further calcula- 
tion and allows the principal solution obtained directly. 


SCANNING CURVES BRIGHTNESS PROFILES 


The relative usefulness scans obtained the three receiving patterns 
without recourse prior knowledge would vary from object object. Some 
appreciation the relative performance may readily obtained reference 
theoretical scanning curves brightness profiles two models small 
source. The models, line source uniform brightness and two equally intense 
point sources, have been chosen since they represent the simplest extreme 
possibilities variation source brightness likely present given 
small region under study the sky, well providing means for showing 
the advantages disadvantages the three antenna patterns. varying 
the extent the models—that is, the length the line object the separation 
the two points—two series scanning curves have been produced for each 
the three antenna patterns. Since the influence the antenna beam 
negligible for uniform sources large angular extent, the series curves has 
been calculated for model angular extents the order the width the main 
exploring lobe the antenna. order make the curves general applica- 
tion for use with antennas different dimensions, unit angle defined the 
cutoff frequency and given the name ‘‘cutoff has been introduced. This 
unit, given the relation and was subdivided into 
intervals .048 for purposes calculation. All the antenna patterns, 
models, and resultant scanning curves have been tabulated and plotted 
these small intervals. Graphical representation suitable curves for illustra- 
tion have been made drawing smooth curves through the plotted points. 
selection the scanning curves obtained with the three antenna patterns for 
various values the model parameters are shown for the line source Fig. 
and for the two-point source Fig. These have been normalized that the 
central lobe always equal unity. The theoretical scanning curves for 
parameter 0.2, with either model, differ slightly from the antenna 
patterns 0), that for the purpose visual inspection either may 
regarded the reference antenna pattern. 

The small differences existing between the various scans are most readily 
revealed measurement the scan features and contain information which 
related the nature the source. Although there are many features the 
curves which might used for description, only the more obvious ones likely 
These are 

(1) the half width the main lobe, 
(2) the visibility the side lobes adjacent the main lobe, 
(3) the visibility the central dip the main lobe. 
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Fic. Calculated profiles line sources various extents obtained with the three antenna 
patterns. 
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Fic. Calculated profiles two point sources various separations obtained with the 
three antenna patterns. 
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graphical representation typical brightness profile scan the 
vicinity the main lobe shown Fig. and arrows indicate the dimensions 
which directly indirectly describe the scan characteristics. will noted 


Fic. Definitions scan profile characteristics. (1) Half width, H.W.; (2) visibility 
first side lobe, (3) visibility central dip, C/B. 


that consideration has been given the zero level and that the strength 
the main lobe taken from the peak the first minimum regardless the 
sign this minimum value. The half width has been taken one-half strength 
the main lobe. Numerical values for the three listed characteristics, half 
width, visibility first side lobe, and visibility central dip, have been 
calculated for the scanning curves presented Figs. and well for 
other curves with intermediate values the model parameter. Figures 
and show these values plotted functions the line width the point 
separation expressed terms the cutoff period. 

Examination these functions shows that the response antenna patterns 
(2) and (3) the angular extent each model increased very similar. 
However, the response the case antenna pattern (1) considerably 
different. the case the two-point model, much more rapid broadening 
takes place from 0.5 1.0, the central dip the main lobe appears 
sooner, and the dip more rapid; the variation the height the first side 
lobe much greater. For the line model similar differences are observed 


TABLE 


Magnitude profile characteristics for various angular extents 


H.W. Central C/B Side lobe A/B 
Point 
Line 
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between the response antenna pattern (1) and the other two, although these 
differences not appear until the angular extent the model greater. 
Table shows the magnitude the three characteristics for values 
0.5, 1.0, and 1.5, derived from the curves Figs. and 
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Fic. Plot the half width the main lobe the profile for various angular extents 
the two models. 
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Fic. Plot visibility the first side lobe the profile for various angular extents the 
two models. 
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Fic. Plot visibility the central dip the profile for various angular extents the 
two models. 


RESOLVING POWER 

The ability given antenna pattern separate two point sources 
expressed most readily terms the appearance the cental dip the main 
lobe, and the case the optical telescope, the least discernible angular 
separation may used measure instrument performance. Many choices 
value for the minimum visibility are possible and have been arbitrarily 
singled one adopting the criterion suggested Rayleigh for solving the 
corresponding problem for the optical telescope. The minimum visibility that 
given two point sources situated that the maximum one single profile 
falls upon the first minimum the other. Antenna pattern (2) will recog- 
nized being identical with the diffraction pattern slit and the arbitrary 
minimum visibility C/B occurs and has value 0.178. For 
antenna patterns (1) and (3) this same value visibility gives the minimum 
angular separation two point sources 0.75 and 0.95 respectively. The 
resolving power reciprocally related this minimum angular separation and 
for patterns (1), (2), and (3) the power 1.33, 1.00, and 1.05 respectively. 
this basis comparison, the sharp-lobed antenna pattern more effective 
separating two sources than the other two. 

However, the presence third source unknown strength and position, 
this antenna with its large side lobes not asset since the interaction 
between the two principal sources may conceal the third one, the third 
source will introduce what amounts error determining the various 
characteristics the scan. For antenna pattern (1), the strength any 
interfering source varies inversely the angular distance the source. 
Antenna patterns (2) and (3), although much broader than (1), not possess 
this degree susceptibility interference since the interference nearby 
source varies these cases inversely the square the angular separation. 


5 15 

— - 
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these two, pattern (3) with its cosine-shaped band pass preferred 
since has the sharper main lobe and the side-lobe level only slightly 
increased. However, appropriate smoothing the observed scan obtained 
with the optimum pattern will accomplish the same result (Barber 


EQUIVALENT WIDTH LINE SOURCE SEPARATION TWO POINTS 


With antenna whose pattern one those considered here, the graphs 
the theoretical curve characteristics may practical use determining 
the nature and size the source being examined. process measuring 
the characteristics observed scan and making comparison with the 
theoretical curves, the most suitable model may chosen. The presence 
absence central dip most useful determining the nature the source 
rather than the half width the main lobe the strength the side lobes. 
Once the nature the source determined, the relative broadening the 
half width the main lobe the profile terms the pattern 
for various angular extents the two models will the most significant, and 
the curves Fig. will special use. For source widths less than 0.5¢, the 
broadening the scanning curves gradual little practical value 
determining the extent the source. 
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Fic. Plot the relative broadening the half width the main lobe the profile for 
various angular extents the two models. 


There are various complicating factors that must considered: sources are 
not often simple the two models chosen for study here; the true antenna 
pattern the observing equipment will have unknown errors; and the presence 
receiver noise will cause errors reading the records. The last may perhaps 
readily allowed for considering the appropriate uncertainty superimposed 
upon the theoretical curves, thus determining the exactness with which any scan 
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feature may read. spite these complications the present development 
for single dimension does show what may expected when sources are 
examined with antennas the same physical extent but different band-pass 
characteristics. 
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APPENDIX 

Methods for obtaining various antenna patterns are suggested operations 
which may performed with equations (1), (2), and (3), while still retaining 
what might called linear relationship between the received power and the 
source temperature; i.e., limiting the investigations essentially the single 
operation addition. 

direct approach the realization this operation involves the construc- 
tion two complete antenna feeders that two distinct antenna patterns are 
formed angles and The powers and delivered the receivers 
are given the equations 


Simultaneous reception these two powers sum, expressed 
the equation 


this equation, new antenna pattern appears formed the sum 
two basic patterns which are separated angle 

The possibilities choice this shift angle are perhaps better revealed 
through repeating the analysis the spatial frequency domain. The individual 
responses the two antennas are given the equations 


where, for example, 

(9a) Ax(s) | A = A(s), 
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Simultaneous reception the two spectra given eqs. (8a) and 
expressed sum the following equation: 


With substitution from eqs. and eq. (10) becomes 
(11) M(s) )A (s) . S(s) 


Thus the effect adding the responses two shifted antenna patterns 
appears the spatial frequency spectrum cosine multiplying factor for 
the original band pass, with the phase term containing the value the shift. 


takes the value one-half the cutoff period, i.e., and 
A(s) isa uniform spectrum, then the cosine term gives the desired cosinusoidal 


band pass the form 


$ 


COs 


Using single antenna, similar result may derived indirectly. The 
response recorded this antenna may used generate second response 
which differs from the appropriate value simply considering the record- 
new variable This follows from eqs. and since shift 
angle appears both sides the equation. Thus the addition two separate 
responses from two antenna patterns replaced the addition two scan- 
ning curves, one derived from the other. practice, this addition could 
readily achieved for fixed scanning rate using time-delay circuit 
derive the delayed scanning curve for its addition the original. 


THE TEMPERATURE-TIME DEPENDENCE THE 
TRIPLE POINT 


ABSTRACT 


The important, but elusive, temperature-time dependency the triple point 
water has been thoroughly investigated triple point cells from two 
sources. During the first days after preparation the cells, the temperature 
was found increase amounts ranging from with the erage 
rise being After the second day the temperature continued rise 
rate about 0.1 per day for about week and finally stabilized. 
practice, ice mantle cell allowed age for about three days before 
the cell used the temperature should reproducible about 

series experiments are described which suggest that this initial temperature 
rise may well due the growth crystals the ice and/or strains the 
freshly prepared ice. The slow rise after the second day could accounted for 
crystal growth. These two possibilities are discussed detail and formula 
relating the temperature crystal size compared with the observed results. 

Tests pyrex cells years old showed that they contain significant 
amount impurities and, therefore, that the segregation impurities during 
the freezing process not likely the cause the initial temperature varia- 
tions. 

the assumption that the above explanations are true, number methods 
eliminating this troublesome initial temperature rise were tested. Since none 
these tests was completely successful, methods extending the usefulness 
old mantles were examined. 

Different methods preparing and using the cells were critically examined; 
the earlier method supercooled freezing was found quite inadequate. The 
effect different thermal bonds the thermometer well and different cell 
environments was investigated. result this work new importance 
attached the standard practice melting the inner layer ice next the 
thermometer well. 

The effect the temperature-time dependency previous measurements 
the difference temperature between the ice and triple points water 
discussed. 


INTRODUCTION 


fixed reference temperature, the temperature equilibrium between 
ice, liquid water, and water vapor (i.e. the triple point) has most important 
role the definition the two main temperature scales. 1954 the Tenth 
General Conference Weights and Measures (Comité International des Poids 
Mesures 1955) adopted absolute thermodynamic (Kelvin) scale defined 
terms only one fixed point, the triple point water, which value 
273.16° was assigned. the International Temperature Scale 1948 
(Comité International des Poids Mesures 1948) recommended that the 
zero the scale realized with the triple point water (+0.0100° 
rather than with the standard ice point (0° C). the range 630° the 
international temperature, computed from the resistance, standard 
platinum resistance thermometer the equation 


received June 29, 1959. 
from the Division Applied Physics, National Research Council, Ottawa, 
anada. 
Issued N.R.C. No. 5384. 


Can. J. Phys. Vol. 37 (1959) 


= 
1230 


BERRY: TRIPLE POINT WATER 1231 


where the resistance the ice point and the constants and are 
determined calibration the steam and sulphur points. Since recom- 
mended that computed from the resistance the triple point 
water, the reproducibility will have bearing every determination 
temperature this range. 

The triple point water was first suggested substitute for the ice 
point 1927 Michels and Coeterier, who found the triple point repro- 
ducible few ten thousandths degree. excellent summary the 
experimental work that followed both the ice and triple points has been 
given Thomas (1941). During this period the ice point was found 
Beattie (1937). the same time the triple point was found 
(1929) and Beattie (1937). While appeared that both points could 
realized with about the same accuracy, the triple point was considered 
susceptible fewer errors because independent both pressure and 
composition the atmosphere. 

The triple point water has been realized glass cells which our cell 
shown Fig. typical example. These cells are filled with gas-free distilled 
water and sealed under vacuum; axial re-entrant well provided for the 
thermometer and usually partially filled with either water mercury 
improve the thermal contact. Details the construction such cells have 
been given elsewhere Barber (1954). 
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Fic. This shows the construction the triple point cells and the four different methods 
ice formation used. All dimensions are given centimeters. 
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Early investigators such Michels and Coeterier (1927), Moser (1929), 
and Thomas (1934) prepared their cells placing them brine bath and 
supercooling the water 10° until spontaneous freezing took place, 
filling the cell with radially elongated needles ice. disadvantage this 
method was pointed out White (1934), who found that after few hours 
the ice tended float leaving interface near the bottom 
the well; this resulted unsteady temperature reading. freezing 
inner mantle ice around the thermometer well, using solid carbon dioxide 
the well refrigerant (Fig. White found could keep cell usable 
for periods week. addition pointed out that this new method 
freezing should concentrate the impurities the water, since they are less 
soluble ice than water, and then when thin layer ice melted next 
the well pure interface would obtained. White found, how- 
ever, that the concentration and diffusion impurities the water was ex- 
ceedingly troublesome and that the inner melting pure ice gave only 
temporary relief. Nevertheless this inner mantle technique was later used 
quite successfully Beattie (1937) and Stimson (1945) and has now become 
standard practice. 

Recent investigations with the triple point cells have confirmed the high 
reproducibility their equilibrium temperature, but have also discovered 
some limitations. Using resistance thermometer, Barber (1954), Stimson 
(1955), McLaren (1957), and the author (Berry 1958) found that the tem- 
perature different triple point cells agreed within Barber and 
Berry showed that pyrex cells and years old still gave the same 
temperature new cells within Because the deuterium oxide content 
water could vary from place place and with the distillation techniques 
used purify the water, Barber felt that the isotopic composition must 
considered order get all triple point cells agree better than 
Both Barber and McLaren found that the triple point temperature some 
cells was low about the day that the mantle was frozen. While 
Barber attributed this effect instrumental limitation also suggested 
that could due strains the freshly prepared ice. Recently Stimson 
(1958) has reported that some his earlier work gave statistical evidence 
that the temperature 3-day-old ice was about higher than 
ice. this laboratory found that each time new mantle was 
frozen cell the temperature usually rose the first 
days* (see Fig. after which remained constant about for periods 
weeks. obtaining this high reproducibility the equilibrium 
temperature also demonstrated that sudden changes due acci- 
dental straining the platinum wire the thermometer could avoided 
careful handling. 

Because the triple point water the most important temperature standard 
and because the initial temperature variations were factor greater than 
the reproducibility the equilibrium temperature, was felt worth 
while try eliminate this initial temperature variation. The first experi- 


*Measurements were usually begun hours after freezing the mantles. 
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Fic. Variations measured two triple point cells with Meyers thermometer 
(from Berry 1958). The first measurement shown with each new mantle was taken about one 
hour after was frozen. 


ments, which were designed accentuate impurity strain effects, led 
understanding the reasons for the initial temperature rise; number 
unsuccessful attempts eliminate the rise then followed. During the course 
this work different methods preparing and using the cells were critically 


EQUIPMENT AND PROCEDURE 

The resistance measurements were made with Leeds and Northrup high- 
precision Mueller bridge using type galvanometer the detector. 
Use effective optical path meters and current reversals allowed 
interpolation between steps the dial for bridge setting 
sensitivity 0.1 ohm per degree was thus possible measure temperatures 
with precision The bridge was given periodic internal cali- 
brations and interpolated bridge corrections were applied the daily measure- 
ments. While this correction compensated for any drift the bridge coils 
relative the internal reference resistor, (i.e. the sum the 
coils) the value was still affected any change itself. The resist- 
ance therefore, was determined periodically terms the ohm 
maintained Canada and found change less than p.p.m.) over 
period years. Since the ohm maintained Canada drifts less than 
p.p.m. year (by comparison with the value presently accepted the 
B.1.P.M. Sévres), the value probably does not change any more 
than few microhms per month. This estimate appears confirmed 
the fact that the previous triple point measurements (Berry 1958) showed 
appreciable drift the equilibrium value over period weeks. 
More details about the bridge and its calibration may found the above 
reference. 

Most resistance measurements consisted observations, current 
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extrapolation from the measurements and ma. Whenever less precision 
could tolerated, for example during the initial temperature rise the cells, 
the above sequence was shortened 1R, 2N, 1N. Before and after 
each measurement the residual resistance the bridge when the dials were 
set zero was determined and its mean value subtracted from the resistance 
readings. All the measurements were made with one Meyers, single layer 
helix, platinum resistance thermometer (S163) which had been annealed 
The thermometer was always placed the triple point cell for 
least minutes before readings were begun. 

Ten cells were used for this work; the age, manufacturer, and difference 
temperature between cells are given Table Unless otherwise specified, 


TABLE 
Description the triple point cells used this work 


Cell Date Feb. 1958, 1959, 
No. Maker manufactured 
N.R.C.* Apr. 1954 
N.R.C.* Apr. 1954 
202 Co.t Feb. 1958 
206 June 1958 +14 


Note: A change in Rrp of 10 wQ is equivalent to 10-4° C, 

*The N.R.C. cells were prepared in this laboratory by Dr. R. S. Turgel. 
+The J. & J. Instrument Co., Silver Springs, Md., U.S.A. 

tA significant amount of air has leaked into cell 25. 


cells were prepared freezing inner ice mantle around the thermometer 
well using solid carbon dioxide (dry ice) and alcohol. The sequence events 
was follows: the first piece dry ice dropped into the practically dry well 
caused dendritic ice crystals form the outside the well the bottom, 
these dendrites then turned into solid glass-like mantle which grew out- 
wards and upwards more alcohol and dry ice were added. The completed 
mantles were tapered, generally being about thick the bottom 
and about 3-5 thick near the surface the water. was found that 
mantles less than 5mm thick the bottom could not used for more 
than day two because hole would form where the ice touched the 
bottom the well. soon the mantle was completed the inner layer 
ice was melted either putting room temperature water the well 
putting chilled water and then warm glass rod the well. The water was 
then left the well increase the heat transfer the thermometer. Before 
each measurement the cell was given sharp rotation check that the 
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mantle was free-turning the well and, therefore, that the inner layer 
water had not frozen. The mantles were maintained good condition for 
periods two months keeping the cells continuously immersed 
bath crushed ice. 
RESULTS 

Test for the Concentration Impurities 

The first line investigation was determine the effect the initial 
temperature rise produced mixing the inner and outer water and melting 
the ice back from the thermometer well. The mixing was accomplished 
inverting the cell several times that the mantle floated and down 
along the well pumped the water into and out the cavity between the 
mantle and well. melt the ice back from the well room temperature water 
was repeatedly poured into the well. 

The results the above tests six different mantles cell 117 are pre- 
sented Fig. the times when mixing melting took place are indicated. 
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Fic. The effect the initial temperature variations cell 117 resulting from mixing 
the inner and outer water and melting the ice back from the thermometer well. The ice was 
melted back about each time this operation indicated. The six tests have been arbi- 
trarily moved along the temperature axis for clarity. 


Throughout this paper the age mantle taken from the time the mantle 
was completed. tests and the mantles were made only three-quarters 
and half their normal length respectively have freer interchange 
between inner and outer water. Since both these short mantles gave normal 
equilibrium values for all the measurements made with them should 
reliable. Some idea the rate temperature change undisturbed cell 
may found the first parts tests and 

each test smooth curve, with about the same rate increase found 
undisturbed cells, could drawn through the points that most them 
fall the curve within their experimental error +2X Even 
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few the sharp changes the curves are real, evident that neither 
mixing the water nor melting the ice back much brought the 
cell immediately into equilibrium. This shows that the inner and outer water 
not contain substantially different amounts impurities (including 
deuterium oxide)* and that the initial rise not due some property the 
innermost layer ice frozen the well. appears, therefore, that the 
initial rise bulk property the ice. 

order absolutely certain that melting the ice back from the well 
and using short mantles did not alter the final equilibrium value the 
tests described Table were made standard mantles that had already 
reached equilibrium. Short mantles were obtained melting those normal 
length down from the surface the water. the results, also presented 
Table II, the only resistance change that significantly greater than the 


TABLE 


Changes the equilibrium value Rrp caused melting the ice mantle away from the well 
and down from the surface the water. All the measurements each test were made the 
same day 


Age of 
mantle 
Test No. days Remarks 
117 still giving correct Rrp spite age. Through usage 
top mantle had melted and ice was back from well 
2mm 
was purposely melted top until only bottom 
mantle was left 
117 Mantle good condition 
Ice was melted back 2mm from well 
117 Mantle good condition 
Ice was melted back 2mm from well. Reading taken 
minutes after previous one 
unchanged. Reading taken hours after first one 
202 Mantle still good condition because had little usage 


Ice melted back 
Inner and outer water mixed 
201 ‘Through usage top mantle melted and ice back from 
well 
was melted back from well 


C); subsequent reading, however, showed that the temperature rise 
was only temporary. evident that even badly deteriorated mantles should 
still give the correct long there continuous sheath ice the 
vicinity the thermometer coil and also that heat conduction the 
thermometer not serious problem. 


Changing the Freezing Rate 

The next line investigation was try freezing inner mantles different 
rates using radically different freezing techniques. Details the four tech- 
niques used are given Table while the graphs the initial temperature 


possibility water undergoing segregation freezing was suggested 
(1958). from data presented Teis and Florensky 
that the initial ice formed might give temperature about low, 


| 
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TABLE III 
Details the freezing procedure for the tests shown Fig. 


Thermal Time 
bond freeze 
well during 
Test Cell freezing Refrigerant minutes 
117 Air Dry ice 
117 Alcohol Liquid nitrogen 
(forced into inner tube) 
117 Alcohol Dry ice 
(seeded top) 
201 Alcohol Brass rod 


cooled liquid nitrogen 
inserted times 


*The thickness given is that at the level of the thermometer coil. 


variations are shown Fig. test the liquid nitrogen was forced down 
the inner two concentric tubes placed the thermometer well. test 
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Fic. The initial temperature variations for cells with mantles frozen varying rates, 
cf. Table for details the freezing The curves are arbitrarily placed along 
the temperature axis. 


the mantle was started the top using piece ice seed and then 
grown down the well. Since all the curves Fig. are similar appears 
that changing the radial rate freezing from 0.2 minute 
eliminating the initial formation dendritic ice near the thermometer 
coil has pronounced effect the initial rise. 

order obtain ice that had been trozen still faster was decided 
try the old method supercooled For this test cells 27, and 201 
were chosen because they had different ages, spread their triple point 
temperatures and were two sources. See Table details. 
The cells were supercooled brine until spontaneous took place 
completely filling them with needles ice. They were then placed ice 
bath and measurements were begun immediately. 
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After hours was found that considerable amount melting had 
taken place and that order preserve the ice, especially the bottom, 
was necessary put the cells back the brine for about ten minutes. This 
resulted the ice needles growing into solid outer mantle about thick 
(Fig. 1C). the loose center ice needles continued melt the outer mantle 
rose until after hours was touching the bottom the well. Finally, about 
twenty hours from the time freezing, was necessary measure the tem- 
perature the outer mantle because the center ice was completely gone 
the vicinity the thermometer coil. With the interface far 
away from the well took the thermometer about one hour reach tem- 
perature equilibrium after being inserted the well and additional 
minutes after each current change. The change resistance due 
measuring current was increased from for inner mantle 
for the outer mantle. 

Figure shows the temperature variation the cells during the first 
days. can seen that all three cells gave practically identical results, 
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Fic. initial temperature variations three different cells with supercooled spon- 
taneous freezes. 


starting out about low and rising equilibrium about one 
day. Most important all was the fact that the equilibrium temperatures 
agreed with each other and with the values for inner mantles better than 
Apparently replacing the cells brine after they were aged hours 
did not appreciably disturb the curves. Although the measurements with 
the outer mantles were extremely tedious make, their values were highly 
reproducible. This pronounced depression the triple point temperature 
during the first day probably explains why Stimson (1955) found super- 
cooled freezes low and why they are generally considered 
less reproducible than inner mantles. 


Test for the Impurities the Cells 

The technique for using outer mantles, once learned, now gave excellent 
check the purity the water; for the mantle were frozen slowly the 
resultant concentration any impurities near the well would alter the 
temperature-time curve and possibly the equilibrium value Rrp. The 
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direction the changes would depend whether ordinary impurities 
deuterium oxide had been concentrated. Employing the same three cells, 
outer mantles were frozen rate 0.2mm per minute placing the 
cells brine bath and using ice seed prevent supercooling. 
After about thickness ice was frozen against the outer wall the cell 
was warmed until the mantle floated and touched the bottom the 
well. The cell was then placed ice bath and its temperature followed 
until equilibrium was reached. further check impurity effects thin 
inner mantle was frozen cell already containing thick outer mantle 
(see Fig. 1D). The layer water between the inner mantle and the well 
would almost certainly mixed with the main body water and therefore 
contain concentrated amount impurities. This latter procedure allowed 
the measurements performed less time and with greater precision. 
The results these two tests well those for the supercooled freezes 
and standard inner mantles are all shown chronological order Fig. 
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O 
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470 INNER MANTLE FROZEN OUTER MANTLE 


TOP TO BOTTOM 4 OUTER AND INNER MANTLE 


FEBRUARY MARCH 


Fic. Variations Rrp measured three different triple point cells with the fourdifferent 
types ice formations shown Fig. The first measurement shown with each new mantle 
was taken about one hour after was frozen. The February measurement cell 117 was 
taken with mantle. 


Allowing for slow drift the measuring equipment, can seen that all 
four types ice formations each cell gave the same equilibrium tempera- 
ture within Also both the outer mantles and the inner and 
outer mantles gave normal temperature-time curves. The absence any 
abnormal effect the tests described above and the fact that old and new 
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cells from both sources have the same triple point within show rather 
conclusively that the total amount impurities the cells extremely 
small. seems certain, therefore, that the initial temperature variations 
our cells, including those for the supercooled freezes, are not the 
concentration and diffusion impurities. The effect can only 
property the ice. 


Fig. most the results obtained with show slow 
linear rise from 0.5 1.5 per day after the comparatively rapid rise 
during the first day two. This most clearly demonstrated the Febru- 
ary mantle cell 201 which gave linear rise from February 
February 25. The resistance drift the bridge and thermometer deduced 
from the values for cell 117 for mantles days old was about 
over the month, not nearly enough account for the observed rise. previous 
work with thermometers the range 100° have detected similar 
long-term drifts from per month. 

the same time was rising visible changes were noticed the 
crystal structure the ice. When mantle was first frozen the crystal boun- 
daries were barely discernible, but after several hours they became quite 
pronounced. Although the crystals were only few square millimeters 
cross section first, they slowly increased size until after week two 
some them had cross sections several square centimeters. That the 
visible crystal boundaries are liquid was demonstrated putting dry ice 
down the thermometer well cell containing old mantle; this caused 
the boundaries disappear and leave clear ice. However, several hours 
later, melting took place, they became visible again revealing the same 
large crystals that had been present the mantle before refreezing. 


Two Explanations for Rising Temperature 

There are two possible explanations for the initial rise terms bulk 
property the ice. The first one, proposed Barber (1954), based 
the annealing strains the freshly prepared ice. Since the Clapeyron 
relation gives per atmosphere the triple point, 
reasonably small strains could cause the observed lowering temperature. 
the case the supercooled freezes the rapid freezing should produce 
greater strains the ice and therefore greater temperature drop. 
unlikely, however, that annealing would take longer than day two. 

The second explanation for the initial rise the triple point temperature, 
suggested Beattie and Stimson (1958), involves the effect crystal size 
and shape the temperature. Since the free energy per mole small crystals 
greater than that the large crystals because their greater surface, the 
vapor pressure small crystals must greater and therefore can seen 
from the Clapeyron relation that the melting temperature must lowered 
the fine state subdivision the solid. The equation relating the 
temperature equilibrium between two phases and the curvature the 
interfacial surface (see, for example, Epstein 1937) pL, where, 


this particular case, 
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the interfacial tension between ice and water ergs 


o 

the heat fusion per gram the triple point, 

where and are the radii curvature the inter- 


face two perpendicular planes centimeters. 
Although the crystals the triple point cells vary considerably size and 
shape shall assume for the purposes calculation that they are spherical 
been reported (Jacobi 1955) but value appears currently accepted 
and used here. The above equation then reduces 


Substituting our observed increase radius approximately 5mm 
order account for additional rise the first day the 
radius would only have change from 0.2 during this time. Con- 
sidering that the actual curvature much greater than that assumed, particu- 
larly the sharp corners the crystals, the effect crystal size and shape 
the right order magnitude account for the results. the case 
supercooled freezes the rapid freezing would produce still smaller crystals and 
therefore cause greater lowering the temperature. Because the relief 
strains and the growth crystals are both accomplished recrystallization 
difficult consider the two effects independently. 

During the course our investigation, out some mantles have 
given temperature-time curves that peaked about long 
days after preparation. These are shown Fig. with the one cell 
also appearing Fig. Both new and old cells exhibited this phenomenon. 
The mantle cell was frozen from the top the well down; the others 
were frozen the standard way. This anomaly not readily explained 
the above theories alone. 
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Fic. four temperature-time curves that showed temperature peaks. The curves 
are arbitrarily placed along the temperature axis. 


/ 
| 


1242 CANADIAN JOURNAL PHYSICS. VOL. 37, 1959 


Attempts Eliminate the Initial Temperature Change 

Based the above explanations number modifications procedure 
were tried attempt eliminate the initial temperature variations. Some 
the results will found Fig. which continuation Fig. 

The following three tests proved quite unsuccessful: 

(a) mantles from the top the well down the hope growing 
larger and less strained crystals near the thermometer coil. See Figs. and 

(b) Freezing the February mantle with cell coupled through oil 
ultrasound generator. was hoped grow less strained crystals. 

(c) the March mantle cell 117 from the top down using 
inner tube hold the dry ice. was hoped distribute the refrigeration 
more evenly and therefore produce less strained ice. 

The next set tests were partially successful that they eliminated the 
large rise during the first day cell 201 but, unfortunately, only reduced 
slightly cell 117. The results are shown Fig. for the following tests: 
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' de 
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MARCH APRIL 


Fic. Variations measured two cells. This continuation Fig. The first 
measurement shown with each new mantle was taken about one hour after was frozen. 
The value Rrp April for the March mantle cell 201 was 25.479536. 


(d) Freezing mantles down the well from ice seed the top 
extremely slow rate putting brine (—1 —2°C) the well. took 
about three hours build 6-mm mantle; even then the very bottom 
remained thin that after the first measurement was necessary 
thicken the bottom dropping small amount dry ice into the well. See 
the March mantle cell 201 and the March mantle cell 117. will 
noticed that thickening the mantle after the first measurement did not 
appear disturb the curve. 

(e) Freezing the mantle first building the top half the mantle with 
dry ice; then freezing about thickness ice the bottom half with 
brine, and finally finishing the mantle with dry ice (see cell 117, April 3). 
This procedure was considerably less tedious than the one test and gave 
the same initial lowering 

Repairing old mantles thickening the very bottom the mantle 


aA. 
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where most the melting takes place. April and April this 
was tried very old mantle cell 201 and left unchanged. However, 
with cell 117 March and April resulted temporary lowering 
respectively. Because the new ice formed had the large 
crystal structure the old ice, rose its value for large crystals more 
rapidly than would have with new mantle. 

Although our observations crystal size during the slow rises 
illustrated Fig. showed direct correlation between size and temperature, 
none the preceding methods freezing new mantle has proved successful 
enough controlling crystal size warrant the additional effort involved 
using them. The graph the March mantle cell 201 well some 
the previous results show that the temperature cell will finally stabilize 
after week two. The difference the behavior cells 117 and 201 during 
the last set tests may due the influence, the strains and crystal 
size, traces impurities possibly the amount etching the glass 
during the initial cleaning the cell. 


Effect not Having Inner Layer Water 

Both Stimson (1955) and White (1934) observed that the temperature 
new inner mantle quite low before the inner layer water formed. 
They assumed that under these conditions the outer interface 
would controlling the temperature and that the lowering was caused 
concentration impurities this outer interface. Since had found 
evidence impurity concentration, even using the outer mantles, this effect 
was investigated again using cells 117 and 201. 

order avoid melting any ice adjacent the well after freezing the 
new mantles, the chilled alcohol was left the well the thermal bond and 
the thermometer was precooled before placing the well. The 
results presented Fig. show that the temperature was low 0.001° 
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within hour after freezing and rose rapidly during the next hour two. 
While the temperature was still low about thin layer ice 
was melted next the well; this caused the temperature jump about 
Asacheck the inner layer water was refrozen again repeatedly 
inserting glass rod precooled until the mantle was longer 
free-turning and the grain boundaries the ice next the well had dis- 
appeared. cell 201 this operation left the measured temperature practically 

more consistent explanation for this large initial rise temperature 
the release internal stress the ice through melting and recrystallization. 
When the inner layer ice was purposely melted there would then 
sudden release the pressure exerted the well the ice. Refreezing the 
inner layer water confirmed previous indications that there difference 
the behavior the two cells. 


Effect Changing the Thermal Bonding Medium and the Cell 

Environment 

order determine whether the thermal bonding medium the well 
has any effect series tests were made with the following media: 
water, alcohol, mercury, and aluminum. The aluminum was the form 
sleeve that fitted closely between the thermometer and the well and was just 
long enough cover the platinum coil. each test the first and last measure- 
ments were made with water the well guard against changes the 
measuring equipment. Unless otherwise specified the height the liquid 
the well was about after the thermometer was inserted. 

the results presented Table the shown the average two 
determinations and should accurate The first time the thermal 


TABLE 


The effect Rrp changing the thermal bonding medium the thermometer well 


Thermal 
conductivity, 


Water 


Pure alcohol +1.4 0.0016 
Alcohol 0.5 water (up hour after mixing) +16 

Water aluminum sleeve* +2.7 2.1(Al) 
Alcohol aluminum sleeve +3.6 

Mercury (about depth) —0.8 0.063 


Water crushed ice down within top 
thermometer coil 


Norte: The values (ARrp) shown in the first column are differences between Rrp with the thermal bond specified 
and Rrp with water. The values (ARima) shown in the second column are the increases in thermometer resistance 
produced by a measuring current of 1 ma. 

*Dimensions of aluminum sleeve are: length 33 mm, outside diameter 11.1 mm, inside diameter 7.4 mm. 


bond was changed from water alcohol, was found high 
one hour after the change. Since addition 0.5 water pure 
alcohol produced similar increase this effect was believed due the 
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exothermic reaction between water and alcohol. This was finally avoided 
flushing the well with alcohol before using thermal bond. 

The results show that water, alcohol, and mercury give the same within 
the experimental error. The increase when aluminum sleeve 
used not likely due exothermic reaction between either 
aluminum and water aluminum and alcohol. From the last entry Table 
appears that few pieces crushed ice should drop into the water 
the well accidentally, would not affected providing, course, that 
the ice not lodged closer than the coil. Although the measuring 
current heat effect for mercury slightly less than for water, water far 
the most convenient use thermal bond. 

1937 Beattie reported that heat leaking downward along the thermometer 
can made negligible immersing the entire triple point cell and part 
the thermometer stem ice bath. also found slight radiation effect 
which was eliminated wrapping the ice bath aluminum 
foil. was decided check for these effects radically changing the cell 
environment from that normally used these measurements, i.e. with the 
cells packed crushed ice stainless steel icebox insulated the outside 
with styrofoam. 

The results the tests are shown Table was found that with the 


TABLE 
The effect Rrp changing the cell environment 


Time Cell environment 


2.10 p.m. icebox. Crushed ice packed top well 
3.00 icebox. Crushed ice only water level cell 
4.04 5215 large silvered glass Dewar 
4.24 671 4.20 cell was taken out ice and left uncovered 
4.35 675 Thermometer was left cell 
4.41 563 4.39 cell was completely wrapped aluminum foil. Still 
4.44 562 out ice 
5.20 Still wrapped aluminum foil but back ice glass Dewar 
5.58 ice box 
Test 
3.22 p.m. ordinary galvanized steel bucket blackened the inside. 
insulation the outside 
4.07 5235 stainless steel Dewar* 
Test 3—Outer mantle cell 
1.50 p.m. Cell completely covered with aluminum foil stainless steel 
Dewar 
2.33 ice box. aluminum foil 


Note: All measurements in each test were made on the same day in the same cell. Unless otherwise specified 
the standard practice of packing the cells in finely crushed ice up to within a few centimeters of the top of the 
extended thermometer well was followed. 

*No difference in Rrp with stainless steel Dewar and icebox was ever observed. 


ice packed only the water level the cell there was still effect from 
heat conduction down the thermometer. Also was not affected measurably 
putting the cell crushed ice different containers such the icebox, 
glass Dewar, stainless steel Dewar, and steel bucket blackened the 
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inside completely wrapping cells containing inner outer mantles 
aluminum foil and packing them ice. However, was affected rather 
seriously the cell was removed from the ice bath. 


Test for Atr the Cell 

this laboratory the criterion for cell containing negligible amount 
air the audibility sharp crack produced the water hammer when 
the cell tipped. With cell this sound longer audible and fact 
large bubbles air can seen the water rapidly inverting the cell. 
Since this cell has temperature only below the others, the above 
criterion seems valid. Also interest the fact that the air cell does 
not significantly alter the curve. 


DISCUSSION 


Two completely experimental determinations the difference tem- 
perature between the ice and triple points water have been made 
Thomas (1934) and Beattie (1937) and two partly experimental determina- 
tions were made Moser (1929) and Stimson (1945); Stimson obtained 
value 0.0100° while the others found 0.0098° With the exception 
Thomas all observers appear have used their triple point cells the 
same day they were prepared. Both Moser and Thomas used supercooled 
freezes; Thomas, however, noticed that the temperature difference was lower 
the first day and therefore used only measurements made least one 
day after the cell was prepared. Stimson appears have used mantles 
the day they were prepared and did not melt the inner layer ice order 
measure the temperature the outside the mantle. appears then that 
most the measurements would affected the time dependency the 
triple point temperature. However, since the measurements would tend 
slightly low, the values finally adopted may more 
accurate than previously thought. 

From the results presented Figs. and can seen that for precise 
intercomparison triple point cells they should prepared the same 
day, allowed age least five days before measurements are begun, and 
comparisons should made only between measurements taken the same 
day. Even then possible that small temperature differences may exist 
between cells due the fact that cells the same age can have different-sized 
crystals. Since this practice was followed only for the 1959 intercomparison 
cells shown Table the 1958 intercomparison shown the same table 
subject greater error. Although our new cells have consistently higher 
triple points than our 5-year-old cells, about not certain 
yet whether the old ones have become slightly contaminated the new ones 
have higher deuterium oxide content. 


CONCLUSIONS 


During the course this work the temperature-time dependence some 
inner mantles was investigated different pyrex cells from two sources. 
was found that the first day after preparation the temperature increased 
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amounts ranging from with the average rise being 
while the second day the temperature increased another 
with the average rise being After the second 
day there was slow rise about per day until finally week 
two later the temperature remained stable within This 
initial temperature rise during the first days thought due either 
strains and/or crystal growth the ice. The long slow rise after the second 
day believed associated with crystal growth. 

series tests, which were designed detect impurities the water, 
were made cells years old. The fact that these tests showed nothing 
unusual and that the old and new cells had triple points that agreed within 
indicates that the total amount impurities our cells extremely 
small. For this reason extremely unlikely that the initial rise due 
the concentration impurities ahead the freezing water. However, 
quite possible that minute traces impurities have some influence the 
temperature modifying the crystal size, their growth rate, the strains. 

the water cell contaminated the technique using inner mantle 
and melting inner layer water cannot give even temporary relief, because 
the time the other factors affecting the temperature have disappeared the 
inner layer water would contaminated. The real value the inner layer 
water having continuous interface close the thermometer 
and the elimination any pressure that the well might exert the ice. 

The old method using supercooled freezes was found result 
temperature that started out about low and rose equilibrium 
about one day. The pronounced lowering the temperature 
sumably due increase strains and/or decrease crystal size, both 
which would occur with rapid freezing. Since the ice was almost completely 
melted after the first day, this method preparation quite inadequate. 

order make mantle that strain free and contains large crystals 
necessary use inconveniently slow freezing rate. One way cir- 
cumventing this problem extend the usefulness old mantle 
refreezing with dry ice and its deteriorated sections; this may 
cause some lowering the temperature but equilibrium reached more 
quickly. Another method use mantle for longer period time. 
have found that they can used for months provided that there continuous 
sheath ice the vicinity the thermometer coil. 

Water was found the most convenient thermal bond use the 
thermometer well and gave the same temperature mercury and alcohol 
within Also was shown that any effect from heat conduction 
down the thermometer from heat radiation was made quite negligible 
completely packing the cell crushed ice. 

order define the temperature the triple point water better 
than consideration would have given not only the deuterium 
oxide content the water but also the crystal size. practice, allowing 
mantle age about three days before using the temperature should 
reproducible Although some the measurements the difference 
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temperature between the ice and triple points water were affected 
this dependence, the value +0.0100° finally adopted 
the International Temperature Scale (1948) probably still within the 
stated limits 0.0002° 
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THE PRESSURE-INDUCED ROTATIONAL ABSORPTION 
SPECTRUM HYDROGEN 


ANALYSIS THE ABSORPTION 


ABSTRACT 


The experimental profiles obtained Part for the pressure-induced absorp- 
tion foreign gas mixtures the region are separated 
semiempirical method into rotational and translational parts. The pressure- 
induced rotational lines are shown have dispersion line shape, modified 
the Boltzmann law. The intensities the rotational parts the spectra are, 
the average, 13% greater than those calculated Part the assumption 
quadrupole interaction alone; the residual experimental intensity probably 
represents the effect overlap interaction. The intensities the translational 
parts the spectra are the same order magnitude the translational 
intensities mixtures the rare gases. Some additional experimental data 
double rotational transitions and the anomalous density variation the 
absorption mixtures are given. 


INTRODUCTION 

Paper this series experimental data the pressure-induced infrared 
rotational spectrum hydrogen were presented and discussion the inte- 
grated absorption intensities was given (Kiss, Gush, and Welsh 1959). The 
experimental results comprised measurements the complete rotational 
spectrum hydrogen the pure gas and mixtures with nitrogen, helium, 
neon, argon, krypton, and xenon; total gas pressures 250 atmospheres 
300° were used, and the measurements were extended when possible 
195° and 85° The spectrum, which lies the region 
showed the expected rotational lines corresponding the selection rule 
+2; general, the lines were very broad with half widths which 
decreased lower temperatures. From values the absorption coefficient, 
a(v), measured across the band, the integrated 
fadv was calculated. This coefficient was shown the 
form Where the density the perturbing gas, except the case 
xenon for which small cubic term, was also necessary. The 
theoretical value the binary coefficient, calculated the basis 
quadrupole interaction collision pairs, was good agreement with the 
experimental value only the case xenon the perturbing gas; all 
other cases the calculated values were appreciably less than those observed. 
When the agreement between the experimental and theoretical values was 
particularly poor, for helium and neon perturbing gases, was noted 
that the intensity ratio the lines did not conform all with that expected 
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from the theory. these cases the rotational spectrum appeared super- 
imposed continuum, the intensity which increased rapidly towards 
longer wavelengths. was suggested that this continuum might collision- 
induced spectrum. This interpretation was confirmed the 
observation similar absorption continua the same wavelength region 
binary mixtures the rare gases (Kiss and Welsh 1959). Translational 
continua are presumably produced the absorption light molecules 
interacting through overlap forces, the energy the absorbed photon appear- 
ing kinetic energy the relative motion the collision partners. 

the further study the intensity the rotational spectrum hydrogen 
necessary separate the rotational spectrum from the superimposed 
continuum. the present communication method effecting 
the separation given, based the fact that has been possible find 
empirically satisfactory analytical expression for the shape individual 
rotation line induced absorption. addition, some further experimental 
data are presented; these concern the intensities double rotational tran- 
sitions and the anomalous density variation the intensity the rotational 
spectrum mixtures. 


THE SHAPE SINGLE ROTATION LINE 
The first step the separation the rotational and translational parts 
the absorption profile obtain analytical expression for the shape 
single rotation line. For this purpose expedient use data obtained 
from low-temperature experiments, since for these the overlapping the 
rotation lines small and the contribution the translational continuum 
less than room temperature. Fig. experimental absorption profiles, 
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Fic. Absorption equilibrium hydrogen 85° the region 300-1000 for 
various gas densities. The rotational Raman frequencies for the low-pressure gas are indicated 
the abscissa axis. 
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vs. frequency, are shown for series densities equilibrium 
hydrogen 85° only the and S(1) lines are present and the over- 
lapping not great. 

The experimental absorption intensities, contain course 
the absorption frequency, which varies greatly over the extent the 
rotational band. the analysis convenient remove this factor and 
transition probability. Fig. the solid line shows this form plot for 
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Fic. Separation the absorption profile for equilibrium hydrogen (85° 118 Amagat 
units density) into two rotational lines the same shape. 


equilibrium hydrogen 85° and 118 Amagat units density. now 
assumed that the two lines, and S(1), have the same shape, separation 
can effected beginning the high-frequency side the band where 
the overlapping least; the result shown the dotted curves Fig. 
The intensity ratio, the two lines 1.06, compared with the 
theoretical value, 1.14, calculated Part the basis purely quadrupole 
interaction. similar separation the lines for normal hydrogen 85° 
gave the ratio, 0.56, compared with the theoretical value, 
0.54. therefore concluded that the assumption similar line shapes 
justified. 

The profile individual rotation line obtained from the above pro- 
cedure asymmetric with the intensity degraded towards higher frequencies; 
the profile has therefore, least qualitatively, the appearance continuous 
distribution summation and difference tones with the rotational frequency. 
this interpretation correct the intensity distribution should the 
Boltzmann relationship, exp(Av/kT), where the intensity 
the high-frequency side the line point from the maximum and 
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the corresponding intensity the low-frequency side. Fig. the plot 
vs. for normal and equilibrium hydrogen 85° gives 


Normal 


Fic. illustrating the Boltzmann relationship, between the 
wings the S(1) line 85° 


acceptable straight line, the slope which corresponds temperature 
The asymmetry the rotational line appears therefore 
expressible the Boltzmann law. 

The most satisfactory analytical expression for the intensity distribution 
the high-frequency wing the line was found the dispersion line 
form, 


where the intensity the frequency vo, the maximum, and the 
value for which the intensity the following, will referred 
the half width the line. evident that the graph 
against should give straight line. The plot for the S(1) line normal 
hydrogen 85° shown Fig. indicates that the dispersion form 
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Fic. Graph illustrating the dispersion line form for the high-frequency wing the S(1) 
line normal hydrogen 85° 
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fact good representation the line shape. The small deviation from the 
straight line Fig. which occurs the extreme part the high-frequency 
wing 250 was observed consistently all the analyses, and 
may indicate real departure from the dispersion shape this region. The 
deviation does not, however, constitute serious limitation the analysis 
procedure discussed below. 

From the above considerations concluded that the shape individual 
rotation line 85° can described fairly accurately combination 
the dispersion line shape and the Boltzmann law. was shown Chisholm, 
MacDonald, Crawford, and Welsh (1952) that the and components 
the split branch the induced vibrational band hydrogen also obey 
Boltzmann relationship. The splitting the branch was therefore ex- 
plained terms the participation the relative kinetic energy the 
colliding molecules the absorption process. view the above analysis 
the wings the pure rotation lines also seem due intervention 
the translational motion the molecules the absorption. There is, 
however, notable difference between the two cases: for the branch the 
most probable change the kinetic energy not zero, whereas for the 
rotation line zero change most probable. This difference between the two 
cases must eventually explicable terms the two different mechanisms 
absorption: overlap interaction the production the branch absorption, 
and quadrupole interaction the case the rotation lines. The explanation 
the dispersion line form for the rotation line must also await theoretical 
treatment the translational energy effect. 


ANALYSIS THE ABSORPTION PROFILES 


The intensity distribution single rotation line 85° has been shown 
above depend three parameters: and the absolute temperature, 
will now assume that the dispersion line form modified the Boltzmann 
law holds for all mixtures and all temperatures with the three parameters 
suitably chosen. The relative values for the different rotation lines will 
taken those given the transition probabilities for quadrupole 
interaction (Kiss al. 1959). Since the temperature known, there are 
effect only two parameters, and the total rotational intensity, adjusted 
fit the experimental data. the example shown Fig. the experimental 
profile for mixture room temperature was accurately reproduced 
variation the two parameters. This particularly favorable case 
since the translational continuum very weak and contributes appreciable 
intensity only the region the line. 

similar analysis was carried out for all the mixtures and temperatures 
for which experimental data were available. Usually only one pressure was 
used since the shape the band changed only very slightly over the pressure 
range. the case and mixtures the separation into com- 
ponents involved some uncertainty since the rotation lines were very broad 
and the translational continuum relatively intense. these cases additional 
criterion was used the analysis, namely, that the difference between the 
synthetic profile the rotation band and the observed profile should 


— | 


CANADIAN JOURNAL PHYSICS. VOL. 37, 1959 


1.0 
Hydrogen-Xenon 
298°K 
Amagat 36.5 Amagat 
02 
200 400 600 800 1000 1200 1400 


Frequency 


Fic. Construction synthetic profile for the enhancement the rotational absorption 
hydrogen xenon: the experimental profile, the component lines the 
synthetic profile, the synthetic profile. 


continuum decreasing monotonically towards higher frequencies. Figure 
shows the final result the separation the observed profiles into rotational 
and translational parts for pure hydrogen and gas mixtures; the 
curves for the mixtures refer the enhancement the absorption the 
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foreign gas, i.e. the intensity due collisions has been subtracted from 
the observed intensities. The ordinate scale used Fig. the specific 
absorption coefficient, giving the absorption intensity per centimeter path 
length per Amagat unit density hydrogen per Amagat unit density 
foreign gas. 

After the division the spectrum into rotational and translational parts, 
possible calculate integrated binary absorption coefficients separately 
for the two parts. Results for the various experimental conditions are sum- 
marized Table For convenient comparison with theory the coefficients 


TABLE 
Experimental and calculated values the absorption coefficients 
Experimental Calculated 

300 0.18 1.19 1.06 155 

195 1.19 1.06 124 


are given the form transition probabilities, and that 
is, with the absorption frequency factor removed, and the units, cycles 
per second per molecule per cubic centimeter hydrogen per molecule per 
cubic centimeter perturbing gas column five theoretical 
values the rotational binary coefficient, calculated Part 
under the assumption quadrupolar interaction alone, are given; the values 
for pure hydrogen differ slightly from those given Part since quantum 
corrections have been taken into account (Van Kranendonk and Kiss 1959). 
The agreement between and now much closer than 
that between and discussed Part The removal 
the translational part the absorption has greatly reduced the outstanding 
discrepancies, for example the case where the translational 
absorption 2.7 times great the rotational absorption. However, the 
values are still most cases greater than and 
probable that the discrepancy (about 13% the average) associated with 
the neglect the overlap interaction the calculation the theoretical 
intensities. This question has been considered detail Van Kranendonk 
and Kiss (1959); they show that, although the rotational intensity mainly 
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due quadrupolar interaction, there small overlap effect arising for the 
most part from interference between the large quadrupolar and the small 
overlap induced moments. The calculation reproduces very well the observed 
passing from normal equilibrium hydrogen 85° using Rosen-type 
wave functions Van Kranendonk and Kiss were able calculate the inter- 
ference effect between the quadrupolar and overlap moments for pure hydro- 
gen. The result was the correct order magnitude but the wrong sign 
explain the difference, this perhaps indication 
the inadequacy Rosen-type wave functions evaluating the angle- 
dependent part the overlap moment. Part the discrepancy between 
and may also caused inaccuracies the calculation 
since some the constants involved, particularly the inter- 
molecular potential constants for mixtures, may error. However, the 
excellent agreement for mixtures 375° and 300° noteworthy. 

The values obtained from the analysis are the same order 
magnitude cm®) those found for mixtures the rare 
gases (Kiss and Welsh 1959). The spectral region which the 
observations were necessarily confined contains only portion the trans- 
lational absorption, which presumably has maximum the lower-frequency 
region. Although quantitative discussion the values 
therefore not very meaningful, some points interest may noted. The 
translational intensity the region 300 decreases rapidly the tem- 
perature lowered might expected. The small value for 
hydrogen must consequence the fact that the overlap interaction 
collisions almost independent the mutual orientation the 
molecules that the induced moment small; the extreme case spher- 
ically svmmetric-like molecules, for example pure rare gas, the induced 
moment zero. the gas series appears have 
maximum for however, the values for and may have 
large experimental error since the translational absorption amounts only 
small fraction the total absorption. 

The last column Table gives the values the dispersion half width, 
the rotation lines, calculated from the slope 1/6) the straight 
line graphs the type Fig. The temperature dependence for 
pure hydrogen and shown Fig. the dispersion half width varies 
and consequently proportional the average relative velocity 
the colliding pair. The total half width the line shows more complicated 
temperature dependence since, from the above analysis, the low-frequency 
wing the line modified the Boltzmann law. The variation the 
dispersion half width given temperature for different collision pairs 
illustrated Fig. 7(6), which plotted against where the 
reduced mass and the Lennard—Jones the pair. The points, 
with the exception those for hydrogen and helium, indicate straight line 
through the origin. Since the factor proportional the average 
relative velocity and measure the distance closest approach, 
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40 


Fic. the dispersion half width, with the temperature, and the 
reduced mass, and the Lennard—Jones ‘‘diameter’’, the collision pair. 


might argued that inversely proportional the time spent 
the region interaction. The linear relationship Fig. indicates that 
inversely proportional the duration the interaction, and can 
construed consequence the Heisenberg uncertainty principle. Since 
the absorption profile shows only small change over the range densities 
used, the density dependence slight. Some measure the density 
effect the half width was obtained from the data for pure hydrogen 85° 
where the rate increase was found about per 100 Amagat 
units density. This the order magnitude ordinary collision 
broadening infrared and microwave spectra nonpolar molecules, and 
adds only small amount the large half width which inherent the 
process induced absorption. 


DOUBLE ROTATIONAL TRANSITIONS THE SPECTRUM HYDROGEN 

85° the and S(1) lines normal hydrogen were observed 
densities 118 Amagat, shown Fig. higher densities, where 
the maxima these lines gave more less complete absorption, two much 
weaker lines were observed the high-frequency wing S(1) (Fig. 8). The 
maxima these weak lines lie close 941 and the calculated 
positions for the double transitions, and which the 
two molecules binary collision perform rotational transitions simul- 
taneously. Since the lines are superimposed the high-frequency wing 
the intense S(1) line, their intensities were difficult measure. particular 
was confirmed that the dispersion line shape for S(1) longer holds 
this region was noted above (Fig. 4); the intensity falls off more rapidly 
than the dispersion line form. Because the difficulty estimating the 
background intensity the binary absorption coefficients for the double transi- 
tions, given Table cannot claim great accuracy; they should perhaps 


considered minimum values. 
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Fic. absorption pure hydrogen high gas densities the region 900-1400 
showing double rotational transitions. 


TABLE 


Integrated absorption coefficients for double rotational transitions 


Experimental 


Theoretical 


The occurrence the double rotational transitions depends the aniso- 
tropy the polarizability the hydrogen molecule. the quadrupole 
moment molecule induces dipole moment molecule the dipole 
moment modulated also the rotational motion molecule the polar- 
izability not completely isotropic. The calculation this effect has been 
given Van Kranendonk and Kiss (1959), using Volkmann’s (1935) value 
the anisotropy the polarizability obtained from the depolarization 
Rayleigh scattering hydrogen. The calculated values are somewhat 
higher than those measured. The discrepancy must partially due the 
inaccuracy the measurement and partially the error 
the value the anisotropy the polarizability. 


THE DENSITY DEPENDENCE THE INTEGRATED ABSORPTION 
COEFFICIENT MIXTURES 
was noted Part that the representation the integrated absorption 


coefficient for mixtures room temperature function density 
was the form, 
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The small cubic term indicates some effect ternary collisions the type 
Since the cubic term was not required the case other mixtures 
the same density range, was proposed that its presence for the 
mixture might due molecular clustering the neighborhood the 
critical point (16° xenon. The molecular clustering would lead 
abnormally high occurrence collisions which some cancel- 
lation the induced dipole would occur, leading negative ternary 
coefficient, 

Experiments the intensity the absorption for mixtures were 
therefore carried out series temperatures 0°, 25°, 75°, and 98° The 
values obtained for the integrated binary and ternary coefficients are given 
Table 


TABLE III 


Integrated absorption coefficients for mixtures 


a1, a2, 


The values are the same within the experimental error for all the 
temperatures. However, the ternary coefficient, has definitely negative 
value and 25°, and too small measured 75° and 98° The 
evidence seems therefore support the suggestion that the relatively large 
cubic term, due molecular clustering the neighborhood 
the critical point. might noted that the isothermals xenon over the 
pressure range used are nearly linear 75° and 100°C, but depart very 
markedly from linearity and 25° 
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THE MAGNETRON GAUGE: COLD-CATHODE 
VACUUM 


REDHEAD 


ABSTRACT 


cold-cathode ionization gauge with axial magnetic field described which 
useful the range above Hg. The gauge operated with anode 
voltage kilovolts and magnetic field 1000 gauss. The relationship between 
ion current and pressure linear the pressure range Hg; 
lower pressures the relationship obtains. the linear region the 
ion current per unit pressure given (amperes) (mm Hg) for nitrogen 
sensitivity about times greater than obtainable with thermionic triode 


INTRODUCTION 


has been shown (Redhead 1958; Hobson and Redhead 1958) that 
Townsend discharge crossed electric and magnetic fields (Penning discharge) 
trapping the discharge space. Moreover, was shown that Penning dis- 
charge gauge could used measure pressures below those measurable with 
gauge (Bayard and Alpert 1950). The Penning discharge 
not limited the X-ray effect which prevents measurements below about 
with gauge. The vacuum gauge described 
these earlier papers was called “inverted-magnetron since the 
anode was located the axis and surrounded cylindrical cathode. The 
inverted-magnetron gauge useful the pressure range Hg; 
its ion current pressure relationship non-linear and given 

where the positive ion current, constant, and the exponent lies 
between 1.10 and 1.15 the normal operating point (anode voltage kilovolts, 
magnetic field 2000 gauss). 

The gauge described herein based similar design principles and has the 
geometry normal magnetron, i.e., the cathode located the axis and 
surrounded cylindrical anode. The magnetron gauge has certain advan- 
tages over the inverted-magnetron gauge for some applications, namely: 

(a) the relationship between ion current and pressure, the normal operat- 
ing point (anode voltage kilovolts, magnetic field 1000 gauss), linear 

(b) the sensitivity the magnetron gauge higher than that the inverted- 
magnetron gauge; 
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(c) lower magnetic field required, thus reducing the magnet size. 
The inverted-magnetron gauge more suitable for measurements below 
Hg. 

The magnetron gauge has sensitivity about times greater than that 
standard gauge (sensitivity (mm electron current 
milliamperes). Since contains hot filament capable measuring 
pressures hydrogen and other active gases which would decomposed 
the hot filament thermionic ionization gauge. 


CONSTRUCTION AND PROCESSING THE GAUGE 
Figure shows schematically the arrangement the electrodes the 
magnetron gauge. The cathode shaped like spool consisting short 


AUXILIARY 
CATHODE 


Fic. Electrode arrangement magnetron gauge. 


cylinder which are welded two circular end disks. The anode short, 
open-ended cylinder which perforated improve gas flow through the 
gauge. Two annular electrodes, the auxiliary cathodes, are inserted between 
the anode ring and the cathode end disks, and are maintained cathode 
potential shield the cathode end disks from high electric fields. The edges 
the auxiliary cathode are rounded and electropolished reduce field emission 
minimum. The gauge enclosed cylindrical glass envelope with 
flattened portion allow insertion between the poles magnet. conducting 
film deposited the inside surface the glass envelope and electrically 
connected the auxiliary cathodes spring contact. This film acts 
electrostatic shield for the cathode end disks. Table gives the salient electrode 
dimensions typical gauge. 


TABLE 
Dimensions magnetron gauge No. 
Dimensions, 

Electrode Symbol 
Anode radius 
Anode length 
Cathode radius 1.45 
Cathode length 


The electrodes are degreased, chemically cleaned, and then fired 1050° 
vacuum for hour. The electrodes are assembled and sealed into the 
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vacuum system rapidly after vacuum firing. The system pumped and baked 
500° for least hours. The gauges receive further outgassing other 
than that caused ion bombardment when the discharge first established. 
After the above treatment there observable evolution gas from the 
gauge the lowest pressures far obtained. The operation and degassing 
the gauge not affected the system let down briefly air and the baking 
process subsequently repeated. 


PERFORMANCE THE GAUGE 


(a) Equipment 
Figure shows the gauge and its associated circuits. The gauge operated 
axial magnetic field about 1000 gauss which provided small 


AUXILIARY ANODE 
CATHODE 


Fic. Schematic diagram magnetron gauge and associated circuits. 


C-shaped permanent magnet. For the measurement pressures greater than 
the anode voltage kilovolts) may derived from simple 
unregulated power supply. For measurements below the anode 
voltage supply must well regulated prevent fluctuations the supply 
voltage coupling the electrometer through the interelectrode capacitance. 
stabilized supply with 0.01% fluctuation for 20% change line voltage has 
been found satisfactory (J. Fluke Co. Model these pressures the 
gauge must shielded from ambient light prevent photocurrents the ion 
collector. The positive ion current measured electrometer the cathode 
lead. The auxiliary cathodes are directly connected ground; thus any field 
emission from these electrodes not measured the electrometer. 


(b) Calibration Method 
continuously pumped systems. The lowest pressures were obtained immers- 
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ing part the system liquid helium (Hobson and Redhead 1958). All 
measurements were made systems constructed Pyrex glass (Corning Type 
7740). All pressures quoted this paper are equivalent nitrogen pressures. 

The magnetron gauge was calibrated directly against Bayard—Alpert 
gauge (Bayard and Alpert 1950) the pressure range where this gauge 
indirect method calibration was employed. 

Figure shows schematically the apparatus used for the calibration the 
magnetron gauge (M.G.) pressures below those directly measurable with the 


) 
VALVE ALVE 
Pa *t U 


Pl 


Fic. Schematic diagram system for calibration very low pressures. 


gauge (B.A.G.). This system consisted two main portions: 
volume which contained the magnetron gauge, B.A.G.-1, and finger which 
could immersed liquid helium; volume which contained B.A.G.-2 and 
was isolated three ultrahigh-vacuum valves (Bills and Allen 1955). Valve 
controlled the conductance between volumes and valve isolated the 
system from the pumps, and valve was used admit gas from reservoir. 

Once this system had been pumped pressure less than Hg, 
valves and were closed. The finger was immersed liquid helium with 
B.A.G.-1 turned off and the M.G. operating. Valve was then opened and gas 
obtained. Valve was then opened slightly and steady flow (Q) gas 
permitted enter volume Then, 


where the conductance valve The steady value obtained was 
determined the input rate gas and the sum the pumping speeds the 
magnetron gauge and the cold finger, thus 


Suat Sr 


where the pumping speed the magnetron gauge, and the effective 
pumping speed the cold finger. Now, 


GAS 
URCE 

{ 
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where the pumping speed the cold surface 4.2° and the 
conductance the tubing between the cold surface and the magnetron gauge. 
has been shown (Hobson 1959) that large compared with these 
measurements, then, 


Gr, 
and since then 


With valve set fixed conductance, proportional Various 
values were obtained adjusting valve and calibration achieved 
plotting the magnetron gauge current against measured B.A.G.-2. 
The value was obtained when was sufficiently high measured 
directly B.A.G.-1. 


(c) Measured Characteristics 

The optimum operating point the magnetron gauge with anode 
voltage kilovolts and magnetic field (B) 1000 gauss. This 
point was chosen because the ion current pressure relation linear and the 
operation most stable under these conditions. Figure shows the relationship 
between ion current and pressure the range Hg, measured 
continuously pumped system direct calibration against B.A.G. These 


CURRENT (AMP) 


PRESSURE (mm Hg) 


Fic. relationship for M.G. No. (with and 
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data were taken kilovolts, 1000 gauss. The measured points 
follow straight line unity slope within +10% for nitrogen, and +20% for 
helium. For nitrogen, 


Uy (amperes) = 9P 


where the pressure (mm Hg). The data for helium were obtained taking 
the sensitivity ratio the B.A.G. for helium and nitrogen 0.209 (Schulz 
1957). The sensitivity ratio the M.G. for the same gases given 


0.15. 


Operation the gauge other values field and voltage gives non-linear 
curves. For linear relation between pressure and current 
must between and kilovolts, and 1000 gauss 10%. 

The operating point must carefully chosen avoid instabilities the 
discharge. Operation away from the optimum operating point results low- 
frequency oscillations over limited pressure ranges and sudden jumps the 
relation. These unstable regions are closely related the 
striking characteristics which will discussed the next section. 

shown Fig. for helium. For pressures above about the 


107 
HELIUM 


Hq) 


Fic. relationship for No. (with kilovolts, and 


pressure relation below this pressure follows power law given 

The data Fig. were obtained the indirect method calibration described 
Section 
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The break the characteristic the magnetron gauge has been observed 
for various voltages and magnetic fields and different electrode dimensions, and 
appears intrinsic this type gauge geometry since the break not 
observed the inverted-magnetron gauge (Hobson and Redhead 1958). 

High-frequency oscillation can observed the discharge almost all 
pressures, the amplitude these oscillations decreases with pressure, and they 
are almost immeasurably small pressures below Hg. The oscilla- 
tions not appear affect the linearity the ion current pressure relation. 
higher pressures care must taken ensure that the oscillations are not 
rectified the electrometer produce spurious readings. 

The striking time for the gauge increases the pressure decreased; 
pressure the striking time about seconds when the gauge 
operated the normal operating point kilovolts, 1000 gauss) 
and when the gauge completely shielded from ambient light. the lowest 
pressure obtained far (about Hg) the striking time may exceed 
minutes when the gauge completely shielded. usually necessary 
irradiate the gauge with ultraviolet light from mercury lamp ensure rapid 
starting these pressures; the radiation removed soon the gauge has 
struck. 

The pumping speed the gauge approximately 0.15 liter/second for 
helium and liter/second for air the normal operating point. One gauge 
constructed entirely titanium showed pumping speed air about 


THE BEHAVIOR THE DISCHARGE 


(a) Electron and Ion Trajectories 

The electron and ion trajectories, the absence space charge, will 
described briefly. Electrons cannot escape from the discharge axial 
direction because the end plates cathode potential. Provided the magnetic 
field exceeds the cutoff field, which the only case interest crossed-field 
discharges, the electrons are trapped radially the magnetic field. Electrons 
will follow hypocycloid path around the cathode while simultaneously 
oscillating back and forth the axial direction. Electrons liberated from the 
cylindrical portion the cathode surface perform one are cycloid and 
return the cathode unless they suffer collision with gas molecule during 
this first arc. The probability electron undergoing collision the first 
cycloidal arc termed cathode escape Electrons 
suffering inelastic collisions are assumed return thermal velocities and 
start their cycloidal motion once more from the point collision. Electrons 
which have escaped from the cathode are liberated within the interelectrode 
space will circle the cathode hypocycloid generated circle diameter 
rolling circle constant radius. The electron trapped indefinitely 
this hypocycloidal orbit until eventually suffers collision with gas 
molecule. each inelastic collision the electron returns near zero energy 
and restarts its cycloidal motion from point closer the anode. Thus 
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electrons can only progress towards the anode successive collisions with 


gas molecules. 
The maximum energy electron its hypocycloidal orbit given 


(1) 


where the electric field 
the anode voltage, 
the anode radius, 

and the cathode radius. 


When less than the ionization potential the gas the electron 
longer capable causing ionization. When 71. Thus for 
ionizing collisions are possible the absence space charge. 

The magnetic field the magnetron gauge considerably below the cutoff 
value for positive ions, thus the ions formed inelastic electron collisions 
move directly the cathode. The transit time for positive ions considerably 
shorter than for electrons, thus the space charge the magnetron gauge 


(b) Striking Characteristics 

approximate theoretical treatment the striking characteristic 
crossed-field discharges has been presented previously (Redhead 1958), 
which assumed that the electric field constant over any one cycloidal 
arc the electron trajectory. The following discussion based 
previously presented theory. 

The conditions for the start the discharge are governed the Townsend 
breakdown criterion which, for non-uniform fields, can written: 


(2) log (14 


where the first Townsend coefficient, the number ions produced 
electron per unit distance advance the direction the electric field, 
and the effective third Townsend coefficient, the number electrons 
re-entering the avalanche process per ion leaving the discharge (i.e., Cyy, 
where the secondary electron emission coefficient the cathode for 
positive ions). 

Figure shows the measured striking characteristics magnetron gauge 
for nitrogen the pressure range Hg. Figure shows the 
corresponding characteristics for helium. 

The diagrams may divided into three regions the lines 
and (see Figs. and 7). the radius which the maximum 
energy the electron its hypocycloidal orbit equal the ionization 
potential the gas, given 


q 
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Fic. Striking characteristics M.G. Fic. Striking characteristics M.G. 
No. for No. for He. 

where the charge-to-mass ratio the electron, 

the electric field the cathode surface, 


and the magnetic field. 


electrons not have sufficient energy cause ionizing collisions. 
region (a), above the (see Fig. 6), electrons are capable enter- 
ing into ionizing collisions with gas molecules throughout the entire inter- 
electrode space. region between the lines and ionizing 
collisions are only possible within the volume bounded the cathode surface 
ionizing collisions cannot occur anywhere the interelectrode space. 

The above analysis ignores the effects space charge and the distortion 
the electric field caused the cathode end plates. can shown that the 
striking characteristics are only slightly altered these two effects. Space 
charge the prebreakdown conditions extremely small and does not build 
any extent until after the Townsend breakdown criterion satisfied. The 
end effects have also been shown small visually observing the start 
discharge observed strike thin lamina, with well-defined boundaries, 
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the middle the cathode cylinder where the electric field greatest and 
where the field distortion the end plates least. The space-charge-free 
potential distribution these gauges has been found from measurements with 
analogue computer which show that the field plane through the middle 
the cathode differs insignificantly from the theoretical field between coaxial 
cylinders infinite length. 

can seen from Figs. and that the line passes through the 
knee the curves and divides the measured striking characteristics into two 
branches. The upper branch the characteristic lies entirely region (a), 
where ionizing collisions can occur anywhere the interelectrode space. This 
branch very closely parabolic, qualitative agreement with the theory 
breakdown crossed fields presented previously (Redhead 1958), which gives 
for the striking voltage the upper branch, 


(4) 
where 
ky (ra—Tx) log (rq rx); 


the striking magnetic field, 
and the effective third Townsend coefficient 
When the dominant secondary processes occur the cylindrical portion 
the cathode surface, given 
(5) P= Cry, 
where the secondary electron emission coefficient the cathode surface 
for positive ions, and the cathode escape probability may approximated 
(see Redhead 1958) 


where the electron mean free path, 


Since nearly constant for the upper branch the striking character- 
istic any changes will caused variations with ion energy; this 
effect quite small. 

The above analysis assumes that the electric field constant over any one 
cycloidal arc the electron trajectory. The maximum radius that 
electron, released from the cathode surface, can reach its first swing outward 
from the cathode has been calculated, and lines constant plotted 
Fig. can seen that the upper branches the striking characteristic lie 
electric field over the cycloidal are very poor approximation for the upper 


branch. This analysis correctly predicts the general shape the upper branch 
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and its movement with pressure, but not capable giving quantitative 
information the secondary processes involved. 

The lower branch the striking characteristics falls region the 
diagram; this branch ionization only possible limited portion 
the interelectrode space surrounding the cathode. The theory gives, for the 
lower branch, 


numerical solution this equation, with assumed constant, shows that 
for gauss, the characteristic rises slowly with However, the 
lower branch will change considerably since approximately constant, 
and varies over wide range. Thus decreases with increasing which 
tends maintain the striking voltage constant. higher pressures, and 
thus larger the numerical analysis indicates that will increase more 
rapidly with found experimentally (see Figs. and 7). 


(c) Running Characteristics 

When the magnetron gauge operated constant anode voltage and the 
magnetic field increased from zero, the ion current observed increase 
rapidly, once the striking magnetic field exceeded, and then become approx- 
imately constant. magnetic field reached the ion current decreases 
instantaneously factor about five and then remains constant this 
new level increased. When the field reduced, upward step ion 
current observed field slightly lower than when increasing The region 
where this step occurs indicated Fig. the region instability. The 
points which the steps occur lie parabola just below the boundary 
between regions and and are only slightly affected pressure. 

Once the discharge has struck, the current limited space charge. The 
effect the space charge, which electronic all pressures below about 
Hg, reduce the field the cathode surface and thus decrease 
the probability escape the secondary electrons. The current through the 
discharge adjusts itself until the field the cathode surface such that the 
effective secondary emission coefficient yXcathode escape probability) 
satisfies the Townsend criterion. the magnetic field increased 
and the operating point moves into region the diagram the space 
charge increases the field the anode above the zero space-charge value 
that electrons are still capable ionizing near the anode surface. some 
value the electric field the cathode becomes small that electrons 
leaving the cylindrical part the cathode never acquire sufficient energy 
ionize before striking the cathode once more. For greater than this critical 
value the discharge maintained only secondary processes from the cathode 
end plates. Since most the ion current strikes the cylindrical part the 
cathode, this critical magnetic field there sharp reduction the number 
secondaries re-entering the discharge and corresponding drop the ion 


current, 


| 
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The magnetic field has little effect the probability secondary electron 
leaving the end plates and entering the discharge since most these electrons 
are moving parallel the magnetic field. Thus becomes independent 
once the critical field has been exceeded. This explains the constancy 
ion current with increasing when B,. The hysteresis caused 
the different space-charge conditions existing either side the region 
instability. 


SUMMARY 


The magnetron gauge capable measuring pressure from 
Hg. has the following advantages over thermionic triode gauge: 

(a) Itis not limited the X-ray effect and thus can measure lower pressures. 

There hot filament cause decomposition hydrogen oxygen, 
evolve carbon monoxide. 

(c) The sensitivity the magnetron gauge the linear region about 
times greater than that standard gauge. 

The pumping speed the magnetron gauge greater than that the triode 
gauge approximately the same ratio the ion currents. Since the range 
operating voltage and field small, the pumping speed cannot readily 
changed the triode gauge, thus the magnetron gauge cannot used 
experiments where the pumping speed must kept small. 

The magnetron gauge most useful the pressure range 
where its sensitivity high and its calibration linear. pressure 
the relation becomes non-linear and the 
gauge used below this pressure must calibrated for the gas used 
the method described Section For measurement these low 
pressures the inverted-magnetron gauge (Hobson and Redhead 1958) 
preferable since there break its calibration curve. 

The measured striking behavior the magnetron gauge qualitative 
agreement with the simple ballistic theory presented previously. 
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FRONTAL PERTURBATION TROPOSPHERIC 
SCATTER PATH! 


ABSTRACT 


2-year study 500 Mc/s radio transmissions over 85-mile path near 
Ottawa has shown that the signal fading rate rises well above the diurnal maxi- 
mum when the transmission path perturbed weather front. The fading 
rate reaches maximum when the upper boundary the frontal zone 3000 
feet above the center the radio path, and the fading rate remains high long 
any part the frontal zone between the surface and 3000 feet. similar 
relationship between frontal position and signal fading rate found for Arctic, 
Maritime, Polar fronts the Ottawa area, but the Arctic fronts provide the 
clearest definition. There some indication that high fading rate also may 
associated with horizontal layers contrasting humidity low levels the 
troposphere. 


INTRODUCTION 

well-known characteristic radio transmission ultrahigh frequencies 
(300 3000 Mc/s) its dependence upon the weather along the radio path. 
When the radio transmitter and receiver are located either within just 
beyond the line sight, the signal transmission loss may estimated with 
reasonable certainty the refractivity profile the local air mass known. 
But the computation transmission loss becomes complicated problem when 
the receiver well beyond the optical horizon the transmitter. There appear 
two reasons for this complexity: First, the long radio path may pass 
through two more different air masses, and each has its own transmission 
characteristics. The second reason even more fundamental, that the 
physical process which the signal travels from the transmitter the 
receiver not clearly understood. least three theories the process 
signal transmission have been suggested (see, for example, Gordon 1955; 
Carroll and Ring 1955; and Northover 1955). The principle that common 
these theories that the radio signal from the transmitter reflected 
scattered within the troposphere and the forward direction, towards the 
distant receiver. 

The results the experiment described this paper suggest new method 
for studying the physical process tropospheric scatter transmission. They 
show that the radio signal path disturbed strongly local weather front, 
whose presence causes the rate fluctuation the received signal increase, 
and the mean amplitude the signal decrease. Further, the perturbation 
the radio path persists long any part the frontal zone pre- 
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ferred height above ground. Thus should possible examine the process 
signal transmission information the physical structure the per- 
turbing front known. 

the purpose this paper describe the experiment that relates the 
frontal position the perturbation the U.H.F. radio path. attempt 
has been made examine the physical structure the transmission path 
detail. The information the frontal positions has been drawn from synoptic 
weather charts and local cloud observations, and the experimental results for 
year have been treated statistically compensate for errors interpolation 
the frontal position. 


RADIO SIGNAL FLUCTUATION 
The experimental radio transmissions were carried out over 85-mile 
path near Ottawa, Canada. 500 Mc/s signal was transmitted from 28-foot 
paraboloidal reflector antenna Chalk River and received 20-foot 
antenna similar construction near Ottawa. Both antennas were directed 
horizontally towards each other. The vertical profile along the transmission 
path, illustrated Fig. shows that the radio horizon within miles 


-1000 
MILES 
TRANSMITTER 
(CHALK RIVER) 
RECEIVER (OTTAWA) 


Fic. Vertical profile along the Chalk River Ottawa path. 


either antenna. The radius curvature the sea level this diagram 
four-thirds the true radius, simulate the refraction standard atmos- 
phere. will seen that the horizon planes intersect about 400 feet above 
the surface near the center the 

Typical scatter transmission occurred between Chalk River and Ottawa, 
that the received signal amplitude fluctuated faded continuously. These 
fluctuations are illustrated the sample records Figs. and The 
amplitude variation was irregular, and the range variation exceeded 
decibels. The rate the signal amplitude generally changed 
orderly way observed equipment designed record the fading rate 
directly and continuously (Day and Psutka 1956). Typical recordings obtained 
with this equipment are shown Fig. The normal variation fading rate 
the 500 Mc/s signal for 2-day interval illustrated Fig. Here, 


details this experiment are given elsewhere published reports (Hay and 
Poaps 1959) and (Hay, Day, and Maynard 1958). 

customary define signal fading rate the number positive crossings the mean 
signal level the envelope the signal, prescribed time interval. this experiment, both 
the prescribed time interval and the interval over which the mean level taken are seconds. 
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(b) HIGH SPEED RECORD SIGNAL WITH HIGH FADING 


RELATIVE 


2000 1600 


(c)RECORD MEAN SIGNAL AMPLITUDE 


Fic. records the time variation the 500 Mc/s signal amplitude 


the maximum fading rate was about per minute, during the early afternoon; 
the minimum fading rate occurred shortly after midnight, when less than 
five fluctuations per minute were 

Soon after the Chalk transmissions began became 
apparent that the normal pattern change fading rate was disturbed for 


brief intervals zero fading rate that are hours apart are due periodic shutdown 
the Chalk River transmitter. 
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FADES PER MIN 


2000 2000 a 1200 
27 MAR. 26 MAR. 1957 
(a) DURING HIGH PRESSURE INTERVAL 


FADES PER MIN 


8 FEB.,1957 7 FEB.,1957 6FEB.,I957 
(b) DURING FRONTAL DISTURBANCE 


Fic. the fading rate the 500 Mc/s signal Ottawa. 


several hours out every days. This disturbance the normal signal 
transmission was shown two ways. The first indication was increase 
signal fading rate, even above the maximum the diurnal change Fig. 
3(a). The record Fig. represents the signal that was received during 
such interval. continuous plot signal fading rate through disturbed 
interval illustrated Fig. Here seen that the fading rate was 
high 100 per minute, and that the fading rate remained high for hours 
longer. 

The second indication abnormality signal transmission illustrated 
Fig. 2(c). This record represents the change mean signal amplitude 
before, during, and after period disturbance. Note that the time scale 
Fig. 2(c) greatly compressed, and that the fluctuations shown Figs. 
2(a) and have been eliminated 2(c) electrical integration over 
several minutes. The important feature Fig. 2(c) the marked depression 
mean amplitude (at the arrow). Such depressions received signal occurred 
whenever the signal fading rate was near maximum during disturbed 
interval. 


DISTURBED TRANSMISSIONS 


the following analysis only those intervals which signal transmissions 
were disturbed have been considered. Further, the criterion selected for dis- 
turbed transmission the abnormally high fading rate, mot the depression 
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mean signal amplitude. Examination the diurnal change fading rate for 
undisturbed transmissions has shown that the maximum fading rate varied 
throughout the year. Table illustrates this variation. Thus, for any interval 
during given month which the fading rate the 500 Mc/s signal exceeds 
the indicated maximum, the transmission assumed disturbed. 


TABLE 
Maximum fading rate for undisturbed transmissions 
Maximum 
Month fading rate Month fading rate 
February August 
March September 
April October 
May November 
June December 


The period from February 25, 1957, February 25, 1958, has been selected 
for the analysis signal fading rate. Although the number useful recording 
days was reduced from 366 306 equipment failures, these occurred 
irregularly, and hence the useful records are representative the full year. 
One hundred and seventy-eight intervals disturbed transmission were 
recorded, and the prominent weather features along the radio path these 
times were examined described below. 

Some indication the relationship between disturbed signal transmissions 
and specific features the weather may found published accounts else- 
where. Prior Hull found that the mean amplitude Mc/s signal 
received distance miles generally was depressed during periods 
precipitation along the radio path (Friend 1945). The rate fluctuation 
the signal increased also the same time. Pickard and Stetson (1947) reported 
later date that prominent signal depression accompanied the appear- 
ance weather front long radio path. The depth depression was 
greatest when the front made large angle with the radio path. Japan, 
has been found that cold weather fronts 75-mile radio path accompanied 
greater variation signal amplitude than did warm fronts (Hirao, Uyesugi, 
and Tao 1952); the extent the variation was greater night than the 
daytime. Experiments radio paths 150 miles long Germany have 
indicated that the depth signal depression varied for the different types 
weather fronts that disturbed the radio paths (Lauter and Klinker 1954; 
Klinker 1956). Randall (1954) the United States has reported signal 
amplitude enhancement prior the depression, when long radio path 
perturbed cold weather front. Thus, increase signal fading rate 
and depression signal amplitude V.H.F. and U.H.F. 
missions have been associated with weather fronts the radio path. However, 
little further information this relationship available. 
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IV. INITIAL ANALYSIS 

Obvious irregularities the local weather were examined first, for their 
relationship disturbances the radio transmission. All the 178 periods 
abnormally high fading rate have been included this analysis. Information 
the local weather was obtained from hourly observations weather stations 
near the radio That precipitation alone has little effect upon the signal 
transmission soon was apparent. Only 12% the maxima fading rate 
occurred when rain fell major part the radio path, and only 
the maxima when snow fell. 

The relationship between disturbed transmissions and air pressure minima 
the radio path was little closer. will recalled that centers low 
pressure generally are associated with strong vertical air currents. Figure 


10 


(a) 


a 


2 40 686 4 2 4 6 8 10 
HOURS AFTER TIME OF PRESSURE MINIMUM * 


INCIDENCE OF HIGH FADING RATE 
n 
a ° 


a2 -6 -4 -2 2 4 6 8 io 12 
HOURS AFTER TIME OF PRESSURE DEPRESSION 


_Fic. Relation between disturbances the 500 Mc/s signal transmission and the drop 
air presssure the radio path. 


compares the time the maximum fading rate with the time minimum 
surface pressure near the center the radio path. all, 102 well-defined 
minima air pressure were recorded during the 306 days signal trans- 
mission. But only these minima were accompanied within hours 
abnormally high fading rate the radio signal. 

higher degree correlation was found between shallow depressions 
surface air pressure and disturbance the signal transmission. These drops 
some millibars air pressure with respect the normal trend lasted 


weather stations, operated the Department Transport, are located Ottawa 
(Uplands), Killaloe (32 miles southwest the center the radio path), and North Bay 
(185 miles northwest the Ottawa site). radiosonde station located Maniwaki, 
Quebec (68 miles north the receiver site). 
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from hours. They were generally accompanied weather frontal 
surface. 306 days, 120 shallow pressure depressions were recorded and the 
signal fading rate reached maximum within hours all but two instances. 
This relationship illustrated Fig. The fact that the peak this 
distribution not the zero time axis may due the separation between 
the radio path and the Killaloe weather station where the air pressure was 
recorded. 

The initial analysis, then, suggests that primarily the structure the 
frontal zone and not the associated precipitation the vertical air move- 
ments the center low pressure that leads disturbance signal 
transmission. This suggestion borne out closer examination the frontal 
heights above the radio path the time signal disturbances. 


FRONTAL HEIGHT AND DISTURBED TRANSMISSIONS 


The continuous changes the height weather fronts above the radio 
path have been deduced from published charts frontal contour 
Information the orderly ascent descent clouds that lay just above the 
front (Grant 1944) has been obtained from hourly observations the weather 
stations near the radio path assist with the interpolation frontal heights 
low levels. convention, the frontal height measured from the upper 
(or warm air) boundary the frontal zone. The full depth the frontal zone 
about 3000 feet, the depth depending upon the type and height the 
front (Godson 1951; Sawyer 1955). Because the frontal zone inclined 
shallow angle the horizontal, the frontal height has been estimated for the 
center the radio path. 

When the frontal data are treated statistically for the whole period the 
experiment, three weather fronts generally are assumed exist above the 
radio path. However, any time this number may decrease according the 
local weather conditions. Polar fronts usually appeared near the surface 
midsummer, Arctic fronts midwinter, and Maritime fronts other times 
the year. Whenever two three fronts occurred simultaneously overhead 
the lower troposphere the mean vertical separation between the upper 
boundaries the fronts was 6000 Thus, when the full interval 306 
days considered statistically more than one type front assumed 
appear within 6000 feet the surface. 

The heights weather fronts above the middle the radio path during 
intervals high signal fading rate are indicated Fig. All the 178 
intervals disturbed transmission are represented these graphs. these, 
154 occurred while weather front was within 6000 feet the surface. The 
charts indicated either front overhead front above 6000 feet 
for the other intervals. should added here that the F.C.A. charts 


7Frontal contour analysis (F.C.A.) charts are issued twice daily the Central Analysis 
Office the Meteorological Branch, Dorval, Quebec. They are based upon the three-front 
model used Canadian meteorologists (Anderson, Boville, and McClennan 1955; Penner 


1955). 
Polar front uppermost, and the Arctic front the lowest when the three appear 


together over the observer. 


we 
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(a) ARCTIC FRONTS 


24 (c) POLAR FRONTS 


INCIDENCE OF HIGH FADING RATE 


° 2 4 6 8 to 12 14 16 18 20 
HEIGHT OF FRONT(THOUSANDS OF FT ABOVE SURFACE) 


Fic. Relation between disturbances the 500 Mc/s signal transmission and the frontal 
height above the center the radio path. 


also indicated weather fronts below 6000 feet during the period the 
experiment when disturbance the signal transmission was recorded. 
Statistics for each the three types fronts found the Ottawa area are 
plotted separately Fig. Fig. 5(c), will seen that the incidence 
disturbed transmissions maximum when the Polar front 4000 
feet. Figure 5(b) indicates two maxima the incidence disturbed trans- 
missions when Maritime fronts alone are considered; one for the frontal 
height 3000 feet and the other for height 10,000 feet. Fig. the 
incidence transmission disturbance reaches well-defined maximum for 
Arctic fronts 3000 feet, and the frequency decreases somewhat lower 
value until the Arctic fronts rise above 6000 feet. For frontal heights above 
this level, disturbed radio transmissions become relatively infrequent. 


VI. REMARKS 


explanation the results given the previous section may based 
upon the known physical structure weather front. Godson (1951) has 
examined the data obtained from many radiosonde ascents Canada, and 
Sawyer (1955) has analyzed the results flight soundings through fronts 
the United Kingdom. They have concluded that there universal model 
front, but that certain features are common most. For example, the 
vertical thickness the frontal zone may between 1000 and 5000 feet 
lower levels the troposphere. Generally, the air the frontal zone 
very dry except levels below 2000 3000 feet. Palmen and Newton (1948) 
have found that the frontal zone possesses strong horizontal wind shear. 
study balloon ascents through well-defined frontal zone has shown, also, 
that the vertical wind shear discontinuous the frontal zone boundaries 
(Reed 1957). 

These features frontal cross section are consistent with the results 
microwave refractometer soundings (Crain, Gerhardt, and Williams 1954). 
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The refractive index profiles show abrupt changes the lapse rate refrac- 
tivity the upper and lower boundaries the frontal zone where the wind 
shear discontinuous. And within the frontal zone where the air relatively 
dry, the lapse rate refractivity approximately uniform but considerably 
different from the lapse rates the air above and below the front. 

suggested here, then, that the upper boundary the frontal zone 
the chief perturbing region the front moves through the radio trans- 
mission path. The remainder the frontal zone presents perturbation 
lesser The frontal zone alters the radio scattering process the 
troposphere that the signal fading rate increased above the normal 
diurnal maximum and the transmission efficiency reduced. Fig. 5(c), 
the Polar front the only one above the radio path when front within 
6000 feet the surface. When the Polar front higher levels, the Maritime 
front and possibly the Arctic front appear below also the radio path. 
Figure 5(b) indicates two preferred heights the Maritime front for 
disturbance the signal transmission; and when this front has risen above 
14,000 feet the incidence disturbed transmissions decreases approxi- 
mately constant value. Thus, appears that the Maritime front the cause 
transmission disturbance rises through the first 4000 feet the atmo- 
sphere, but the disturbance when the Maritime front 10,000 feet due 
the Arctic front below (now within 4000 feet the surface). Finally, 
Fig. suggested that the high fading rate results from perturbation 
the radio path the Arctic front Perturbation the radio path 
becomes minor when the Arctic front higher than 6000 feet. 

Figure illustrates the suggested explanation radio path disturbance. 
consistent with the experimental results given above assume that only 


Fic. section through the radio path and the perturbing front. 


reasonable expect more abrupt change the lapse rate refractivity the 
upper boundary the frontal zone than the lower boundary, because the air warm (and 
relatively moist) above this zone but cool (and drier) below. 

Maritime front now 6000 feet above the Arctic front, and the Polar front 6000 


TMTR RCVR TMTR RCVR 
(a) POSITION FOR MAX, PERTURBATION (b) HIGHEST PERTURBING POSITION 
feet above the Maritime front. 
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the lowest the three fronts causes significant The diagram 
represents vertical section through the radio path when occupied 
weather front. The front, which extends perpendicularly the plane 
the diagram, somewhat idealized here for the sake illustration, and 
only the lowest front considered. The vertical depth the frontal zone 
assumed 3000 Figure represents the position the front 
when the signal fading rate maximum. Here, the upper boundary 
about 3000 feet above the surface. Because the fading rate remains high 
level for frontal heights between the surface and 6000 feet, only for frontal 
positions above that shown Fig. that perturbation the radio path 
becomes negligible. Note that the lower boundary 3000 feet above 
the surface. The frontal motion the right about miles per hour 
the lowest 3000 feet the troposphere, and especially the 3000-foot level, 

appear critically important radio transmissions over 85-mile path. 
This description frontal perturbation has assumed simplified picture 
the troposphere. That other sources perturbation may present 
suggested the fact that the incidence transmission disturbance does not 
vanish Fig. when all fronts are higher than 6000 feet. Twenty-four 
intervals this type were recorded during the experiment. examination 
the radiosonde data from Maniwaki for these intervals showed prominent 
irregularity the lapse rate temperature and relative humidity within 
5000 feet the surface instances. These irregularities generally represent 
horizontal layers contrasting humidity the air. Since horizontal stratifi- 
cation this type often follows the curvature the earth’s surface for 100 
200 miles (Bauer, Meyer, Prohaska, and Wilson 1956; Crain 1955), 
probable that such layers also contributed the disturbance signal trans- 


missions. 
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LIMITING PRECISION OPTICAL 


ABSTRACT 


The fundamental limit the precision setting fringe peak due 
photon noise evaluated for certain class interferometric methods terms 
parameters characterizing the source, the interferometer, and the detector. 
shown that changes path difference should detectable with 
observing time second, using only modest equipment. Some the experi- 
mental conditions required attain this precision are discussed. 


INTRODUCTION 

The corpuscular nature radiant energy imposes limit the precision* 
any optical measurement that occupies finite observing time. Thus, even 
ideal experiment which all the available photons are effectively used 
there will statistical spread the results successive 
measurements the same quantity, that is, the precision finite. possible 
calculate the maximum theoretical precision expected any particular 
method measurement, and comparison this with the observed precision 
will indicate how successful the experimenter has been eliminating least 
the fluctuating part possible perturbing effects that can influence the 
results. The steady part, bias, such effects reduces the accuracy* but not 
the precision the measurement. principle these perturbations can 
eliminated, leaving only the fundamental quantum limitation precision. 

Some the most careful interferometric measurements recent years 
have been made connection with international program that has cul- 
minated resolution, presented the International Conference 
Weights and Measures, redefining the Metre terms the vacuum wave- 
length the 6056 line The wavelength this radiation has been 
determined terms the cadmium red line number national labora- 
tories. The over-all accuracy this comparison, judged the agreement 
the results obtained the various laboratories, about The precision 
reported for some the individual results was somewhat higher. 

The question arises how close have come the theoretical limit 
precision measurements this type. partial answer this question, 
certain class interferometric methods using modern techniques de- 
scribed, and analyzed detail Part I]. The maximum precision that can 
expected evaluated terms parameters characterizing the source, the 
interferometer, and the detector. number conclusions are drawn from the 
functional form the expression, and numerical examples are given. Since the 
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*In this discussion, precision defined the reciprocal the relative standard deviation 
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theoretical precision turns out much higher than the accuracy obtained 
present, Part III devoted consideration some the experimental 
conditions that must established attain this precision. 

The general form the apparatus considered sketched Fig. 


Fic. Schematic diagram experimental arrangement. 


the optic axis. Light from each point passes through the interferometer 
parallel beam and emerges from with diameter The interferometer, 
which may Fabry—Perot etalon, two-beam-type like the Michelson, 
has its reflecting surfaces adjusted normal the axis and therefore 
produces fringes equal inclination centered The interferometer also 
has arrangement for making small adjustments the order inter- 
ference, making small changes the optical path difference. 

the normal optical path difference the interferometer the order 
interference the center the ring system has the value 


(1) 2to 


where the wave number the center the spectral line. (It assumed 
that the line symmetric.) changed, the intensity the light passed 


Fic. Variation flux with order interference when small aperture used. 


has maximum, that is, set the fringe peak. The fractional part 
follows immediately from the change necessary make the setting; 
The integral part can obtained several methods and will 
not considered here. 

The practical method setting the fringe peak compare the in- 
tensities order numbers which are located symmetrically 
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about and separated fixed interval. The criterion for making the setting 
seen the detailed analysis below that, practice, the aperture should 
not very small, but has optimum size. According our earlier dis- 
cussion, the required intensities cannot measured exactly, and thus there 
will finite minimum departure from the condition that can 
detected above the unavoidable noise level the detector. The maximum 

instrument using essentially the above technique with 
etalon was described Baird (1954). The setting was made visual detec- 
tion minimum flicker with c.p.s. scan between and obtained 
varying the air pressure with motor-driven piston. later version with 
photoelectric detection was used measure the shift some mercury lines 
with carrier gas pressure (Baird and Smith 1957). Terrien has applied 
the method Michelson interferometer, making static readings photo- 
multiplier current first and then has calculated the precision 
expected his instrument when the limiting sources noise were (a) 
shot effect the photocathode current and lamp fluctuations. the 
analysis which follows considered that only the former fundamental 
limitation, the lamp fluctuations being principle eliminable. this basis, 
the limiting precision will evaluated for both two-beam and multiple-beam 
interferometers with plane reflecting surfaces, and obtained form that 
displays its dependence the parameters the apparatus. 


ANALYSIS 
Each point the plane (Fig. characterized its angular dis- 
placement from the optic axis seen from the center angle 


the order interference will less than that the center amount 
where 


(1) and (2) obtain immediately 


has focal length the area the lamp used illuminate orders 


(3) 
The solid angle over which light collected from each point the source 
(4) 


the source has radiance (watts/metre? steradian) the light given 
line, the power incident the interferometer which available produce 
the interference pattern for orders is, from (3) and (4), 
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Since this linear the power available for each order the same, and the 
power emerging from the interferometer each order 


where the transmittance the instrument averaged over one order. 
For the Michelson interferometer 


where are the transmittance and reflectance the beam splitter, and 
the reflectance the mirrors. 

practical instrument will somewhat less than its maximum 
theoretical value 0.5, but the assumption this value will not significantly 
overestimate the performance. For the Fabry—Perot etalon 


where are the transmittance and reflectance each etalon plate. Because 
the multiple reflections, the effect absorption the reflecting films 
more important than the Michelson, and the actual values must 
used the above expression obtain The power (equation 
distributed over the order pF(n)dn where normalized such that 


mot} 
F(n)dn 


make setting the fringe peak observing time the fluxes re- 
ceived from each side the fringe during are compared. Since the fluctua- 
tions these two signals are independent, the root-mean-square fluctuation 
their difference will the same that the sum, i.e. the noise will 
determined the total number photons received during time 
given 


where the energy one photon and are the order numbers 
the edges the band passed the aperture (Fig. 1). (We need only 
consider the signal from one side the fringe since the signals from the two 
sides are equal when the setting has been made.) With only small error 
may assumed that the source steady, photons are emitted randomly, 
and thus the r.m.s. fluctuation will the presence this 
noise take the smallest detectable departure the setting from the peak, 
that value which produces difference signal just equal 
the r.m.s. fluctuation where 
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the photons that fall the detector, only fraction produce count- 
able output. has been called the quantum yield (Rose 1948) and the detective 
quantum efficiency (Jones 1958). will equal the quantum efficiency 
the primary photon detection mechanism only there are appreciable 
noise sources between the detector and the final indicator.) calculate 
must now equate Using equations (1), (5), (6), and 
(7), and assuming spectral line width the maximum precision with 
which can determined 


2 1/2 


and the fraction order occupied the line width Ag, i.e. 

The next problem find values and that make maximum. 
Gaussian line shape will assumed, but the general results will not 
sensitive line shape. For the Michelson interferometer 


where the visibility the fringes, given 


For the etalon 
k=1 


This expression was derived Krebs and Sauer (1953). 

Numerical investigation several representative curves revealed that 
was largest when the interval a(= embraced most the side the 
fringe, and that value within 20% the maximum could 
obtained all cases choosing and such that 

These assumptions greatly simplify the calculation and retain enough 
accuracy for our purpose. follows immediately from (9) that 


§ 
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Then 


For the Michelson interferometer, simple calculation shows that 0.205. 


Taking 1/2, becomes 


0.565. 

The calculations for the Fabry—Perot etalon could not done readily 
because the complication the function Numerical values were 
determined from the tabulated curves given Krebs and Sauer (1953), and 
are shown Fig. function for several values for dielectric and 


.02 04 O7 2 4 
Fic. Variation with dielectric films, aluminum films; 
Michelson: 


aluminum reflecting films. The individual curves for aluminum are the same 
shape those for dielectric films with the same value but the heights 
are reduced because the higher absorption. This effect particularly severe 
for large The circles these curves indicate the point which the line 
width equal the instrumental width. The curve for the Michelson interfero- 
meter also shown the same figure. 


02 04 .O7 4 7 | 

73% 

61% 


HANES: OPTICAL INTERFEROMETER 1289 


seen from Fig. that the maximum attainable precision about the 
same for the Michelson interferometer and the unless 
necessary use aluminum films. (Aluminum has very little dispersion 
phase change reflection, property that useful when comparing wave- 
lengths different parts the spectrum.) this case, the 
suffers more because the multiple reflections. 

The path difference etalon length which the optimum precision 
obtained relatively short, ranging from about 1/14 1/7 depend- 
ing the instrument. The smallest obtained from conventional sources 
given the lamp operated the triple point nitrogen. For this 
source ranges from about cm. The optimum reflectance used 
the Fabry—Perot surprisingly low, about 70% for dielectric and 55% for 
aluminum films. 

important fact displayed equation (8) that the precision attainable 
depends the characteristics the spectral line only (There 
indirect effect arising from the wavelength dependence for the detector 
used.) For given line shape, directly related the power emitted 
per unit spectral interval the center the line. Thus any arrangement for 
producing narrower line, such cooling the source, the use atomic 
beam, will lead loss attainable precision unless the peak intensity 
maintained. practice, smaller value highly desirable, since the 
interferometer used larger path difference where both and are larger 
for the same precision. The larger allows relaxation the tolerances 
steadiness the source, uniformity illumination, flatness, and parallelism 
the interferometer mirrors and, course, larger fringe fractions are easier 
measure. These practical gains may lead acceptance source with 
smaller limiting precision (i.e. its smaller half width allows higher 
precision actual experimental practice. should noted that accepting 
loss only factor the theoretical maximum precision allows in- 
crease approximately 0.6. The same practical advantages accompany 
this change those resulting from the use narrower line. 

example, the value equation (8) will calculated for the 6056 
line discharge lamp operated the triple point nitrogen. MKS 
units are used throughout. For this line (Engelhard 1958), 
photomultiplier with trialkali cathode, such the 
series EMI, has this wavelength. Using etalon with 
5-cm spacer and reflecting plates diameter coated with dielectric films 


For the green line Megger’s type lamp, the value about 
five times larger with the same photomultiplier. The maximum precision for 
this line therefore about The value the precision obtainable 
with the krypton line confirms the estimate made Terrien for his 
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Michelson interferometer. His calculated precision was 10° for 
cathode current about 10° electrons per second. The value used here (see 
below, Section electrons per second, which gives precision about 
times greater, i.e. about 3X10". 

experimental confirmation (8) afforded some figures kindly 
supplied Smith. These results were obtained while measuring the 
shift the krypton orange line (Baird and Smith 1959) using later version 
the apparatus described Baird and Smith (1957). Putting the con- 
stants relating his apparatus, equation (8) gives terms 
detectable wavelength shift this The estimated precision 
single setting was about Since shift just equal the r.m.s. fluctuations 
would not reliably detected these values are considered good 


agreement. 


EXPERIMENTAL CONDITIONS 


this section some the experimental conditions that must established 
allow approach the above theoretical precision are evaluated. Num- 
bers used for illustration relate the first example described above using the 
krypton orange line. The material this part intended give some idea 
the problems involved attaining such high precision. should stressed 
that the attainment accuracy measurement the same order 
vastly more difficult task, since there are very many sources systematic, 
but usually constant, errors that are not even mentioned here. 


(a) The Fringe Fraction 

The analysis Part shows that change fringe fraction from 
the setting the peak produces noticeable output from the detector. 
(It assumed, course, that sufficient amplification provided between 
the detector and the final indicator show the detector noise.) the example, 
the mechanism for adjusting the fringe fraction must fine enough that 
changes this size can made and read. This requirement not prohibitive 
since many experiments require only small range adjustment 


(b) Stability the Interferometer 

attain the precision the example, changes other than those deliber- 
ately introduced make the setting must less than parts (i.e. 
1/40 cm) over any period time taken for measurements. This re- 
quires high order mechanical stability, and the temperature must 
held constant about even with invar fused quartz construction. 


(c) Lamp Stability 

straightforward method setting the fringe peak measure the 
flux from one side the fringe for time 7/2, shift the other side the 
fringe for and compare the results. evident that the difference 
average flux leaving the source two successive intervals must less 
than the variation the number counted events that are expected from 
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quantum fluctuations. The relative r.m.s. fluctuation output due quantum 

The assumptions (9) and (10) led value the integral (6) 0.205 
for the Michelson. The same value, within 5%, was found all the cases 
considered for the Fabry—Perot well. Thus may calculated from 


the example above, 1.4X10° and Thus the 
stability the light output from the source, described above, must 
better than 0.01%. 

might mentioned this point that the photocathode current this 
example, amperes, much larger than the dark 
current amperes). Thus the shot noise the dark current neg- 
ligible and, ignoring the small amount noise introduced the electron 
multiplier, the quantum efficiency the photocathode may taken for the 
detective quantum efficiency the tube. 

Another method setting the fringe has been used the National 
Research Council (Baird and Smith 1957). this method the order changed 
relatively rapidly (15 cycles per second) between and sinusoidal 
scan, the mean value being adjusted give null the output 
synchronous detector operated the scanning frequency. The lamp fluctua- 
tions that affect the setting precision are narrow frequency band the 
scanning frequency. This method will then less sensitive lamp fluctua- 
tions the spectrum these fluctuations has decreasing quadratic content 
higher frequencies. 

brief experiment investigate this point was performed. Light from the 
source being studied fell the cathode photomultiplier tube. The current 
from the anode passed through load resistor and microammeter. The alter- 
nating voltage across the load was amplified, passed through narrow band 
filter tuned and read r.m.s. meter. From previous measure- 
ments this phototube, the r.m.s. meter reading function average 
anode current was known for steady source. This course the amplified 
shot noise the photocathode current. the lamp fluctuations are large 
enough, higher meter reading will expected. Since the lamp fluctuations 
can considered modulation the incoming light, the relative r.m.s. 
value will independent the photomultiplier current, whereas the relative 
shot noise will vary inversely the square root the current. using the 
whole spectrum from the source, much larger currents than are normally en- 
countered could obtained, and therefore lamp fluctuations that would 
normally buried the shot noise could measured. 

For the lamp operated the example above from unregulated 
power supply, the relative r.m.s. fluctuation c.p.s. pass band c.p.s. 
was Thus the fluctuations output this lamp will not appreciably 
reduce the precision when the rapid scanning technique used. With the 
static method setting described first, the important fluctuations are those 
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occurring successive 1/2-second intervals. These are very much larger, 
approaching some cases, although the exact r.m.s. value was not deter- 
mined. Considerable stabilization the source would therefore necessary 
this case. Similar measurement made electrodeless lamp excited 
from unregulated oscillator about 300 Mc/s showed relative r.m.s. 


(d) Illumination 

Lack perfect flatness and parallelism the reflecting surfaces the 
interferometer will cause variations over the field. The effective value 
will then depend the illumination the field. 

The center the aperture (A, Fig. corresponds and (Fig. 2), 
and since these are not symmetrically located relative the fringe peak 
ny, illumination over the aperture will cause apparent shift 
the peak. Approximate calculations these two effects show that variation 
the order the illumination may lead errors few parts 
certain types measurements. However, most experiments the 
relative illumination will constant, and loss precision will occur. 


IV. CONCLUSIONS 

The principal results this analysis are given equation (8) and Fig. 
They show that the precision obtainable with the 
Michelson interferometer are the same order, and confirm 
calculation for Michelson interferometer with photoelectric detection 
(Terrien shown that contrast the requirements for high 
resolution, optimum precision setting fringe obtained with relatively 
short path difference, relatively large aperture the plane the Haidinger 
rings, and, for the mirrors quite low reflectance. According 
the theory, the only important property the source the radiance 
the light given line divided the line width. However, practice 
small line width gives important experimental advantages. additional 
consideration practice the wavelength dependence the detective 
efficiency available detectors. 
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THE SPECTRUM Li* 


HERZBERG AND Moore? 


ABSTRACT 


The spectrum the Lit ion has been reinvestigated under high resolution 
the region 10,000 1000 with view determining the Lamb shifts the 
ground state and first-excited state. The hyperfine structure the line 5485 
has been studied under improved resolution both with the natural isotopic 
mixture and with 96% Li® sample. For many the higher levels the hyperfine 
structure larger than the fine structure and for large and even larger 
than the separation. The coupling conditions these cases are 
discussed some detail and theoretical intensity distributions the combined 
structure and hyperfine structure are presented for comparison with the 
partially resolved observed structures. The selection rule for the quantum 

the basis this detailed analysis table energy levels Li* presented 
including, where possible, the hyperfine splittings. The absolute accuracy the 
energy values only because the error introduced the resonance 
lines below 200 not remeasured the present work. However, the 
relative accuracy the levels with reference the level considered 
within +0.05 most cases. 

The limits the observed series have been determined with accuracy 
and lead ionization potential 


610079., 


with error caused entirely the error the resonance lines. The 
term values (ionization potentials) the 28S and states 


134044.19+0.10 

120008.30+0.10 
are much more accurate. The the observed ionization potential 
and that calculated from relativistic quantum mechanics —8.0+3 which 


agrees within the rather large limit error with the predicted Lamb shift 


INTRODUCTION 


Most the experimental and theoretical work the Lamb shift refers 
the one-electron systems However, the problem Lamb shifts 
two-electron systems has been attacked recently. Kabir and Salpeter (1957) 
and Sucher (1958) have carried out extensive calculations leading pre- 
dicted value for the shift the ground state and estimates for 
and higher ions this series. Experimentally, the Lamb shift can deter- 
mined the difference the observed ionization potential and the ionization 
potential calculated the basis the Dirac theory (neglecting quantum- 
electrodynamic effects). Both quantities have now been obtained with con- 
siderable precision the case (Herzberg 1958; Pekeris 1958) 
difference agrees very satisfactory way with the predicted Lamb shift. 
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Since the accuracy neither the predicted nor the observed values for the 
Lamb shift very high, appears desirable obtain experimental 
Lamb shifts few other cases energy levels two-electron systems. For 
this reason, work the spectrum Li* was initiated this laboratory. 

experimental determination the Lamb shift possible not only for 
the ground state but also for the metastable state two-electron systems. 
theoretical (Dirac) value for the energy the level has been 
calculated Pekeris (1959) and that for being calculated Hart 
and Herzberg (1959). For precise experimental value the energy 
the level was implied the earlier work the ground state He, namely 


Previous work the spectrum has been summarized Moore 
(1949). Term values were obtainable that time only within and the 
ionization potential was estimated have error +25 Such low 
accuracy clearly insufficient determine Lamb shifts and for this reason 
the present measurements were undertaken. has been possible obtain 
the term values large number excited states (including with 
accuracy +0.10 but definitive term value for the ground state (i.e. 
for the ionization potential) has not yet been obtained, although considerable 
improvement over the previous accuracy has been accomplished. For precise 
determination the ionization potential the resonance lines 199.3, 178.0, 
will have remeasured with much increased accuracy. Unfortu- 
nately these lines are spectral region that have not yet reached with 


our present equipment. 


LIGHT SOURCE 


Because the high ionization and excitation potential Lit (the highest 
any singly ionized system) its spectrum not easy obtain, particularly 
since the low vapor pressure and high chemical activity add further 
difficulties. This explains the scarcity previous work the spectrum 
and the low accuracy obtained. 

Three different kinds sources were tried the present experiments, (a) 
high-frequency, high-voltage discharge with external electrodes, sliding 
spark (see Bockasten 1955), and (c) hollow cathode tube. The high-frequency 
discharge did not yield any lines. The sliding spark gave good intensity 
but apart from other difficulties the lines were quite broad, presumably 
account Stark effects due interionic fields. The hollow cathode was 
far the best source both for its intensity and for the narrowness the 
lines produced and was therefore chosen for the present investigation. The 
original design Schiiler (1926) was slightly modified. The cylindrical cathode 
length and 1.4-cm diameter was made molybdenum, which was 
found withstand the chemical activity lithium high temperatures 
fairly well. The internal diameter the cathode was 1.1 cm. One end was 
closed while the other had opening 0.6-cm diameter. These dimensions 
approximate those recommended Risberg (1956) his work and 
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The molybdenum cylinder was surrounded close-fitting copper 
jacket external diameter This jacket ensured more uniform 
distribution temperature along the cathode than one would otherwise 
obtain. The cathode with its axis horizontal was supported vertical 
thin-walled steel tube which insulated the cathode thermally and the same 
time, when necessary, could used cool the cathode stream cold 
air. front the opening the cathode was small iron ring which served 
the anode. Although the same (ground) potential the whole metal 
vacuum chamber, this anode appeared stabilize the discharge and lead 
increased intensity the lines. 

For the operation the discharge, helium pressure was used 
carrier gas. d-c. potential 350 volts was applied and currents 
300 were used. charge about 0.3 lithium was introduced into 
the cathode and was sufficient for exposures hours. The cathode 
could heated 1000° directly the discharge. separate furnace, 
used Schiiler, was not found necessary evaporate the lithium. 

The superiority the hollow cathode source over the sliding spark source 
shown Fig. which reproduces spectra obtained with the two sources 
the fine structure the line 5485 


OBSERVED SPECTRA AND STANDARDS 


The spectra the photographic infrared, visible, and near ultraviolet 
regions (8600-2000 were taken the first second orders the N.R.C. 
21-ft and 35-ft grating spectrographs. The spectra the vacuum region 
(2000-1180 were taken the fourth and fifth orders the 3-meter vacuum 
spectrograph. 

Above 2000 and lines from separate hollow cathode tube, photo- 
graphed the same order the lines, were used standards. order 
eliminate any possibility shifts between the Lit lines and the comparison 
spectrum (on account changes temperature barometric pressure), the 
latter was exposed partly before and partly after the main exposure. Figures 
and show sections these spectra. 

the vacuum ultraviolet, order obtain comparable wave number 
accuracy, imperative obtain higher wavelength accuracy than the 
visible and near ultraviolet regions. Therefore, internal standards occurring 
the same source the spectrum measured and the same wave- 
length region are desirable. Such internal standards near the five lines 
the vacuum region are lines and derived from the com- 
bination principle Edlén and Risberg (1956) and Risberg (1955) and the 
resonance lines near 1200 measured Herzberg (1958). The former 
were produced the source adding small amounts calcium mag- 
nesium the hollow cathode, while the latter were obtained mixing 
trace nitrogen with the carrier gas. example, Fig. shows the 
lines 1493.0 and 1653.1 together with two Cat lines. Other Ca* lines 
the second spectrogram, well lines the first, are beyond the 
region reproduced. The spectra shown were taken the fourth order the 
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vacuum spectrograph and were separated from the overlapping second and 
third orders the use foregrating described Douglas and Herzberg 
(1957). spectrum from the hollow cathode tube mentioned super- 
imposed second order (by suitable rotation the foregrating) this 
serves only for the determination the shape the correction curve whose 
absolute position determined the and lines. 

The exposure times varied considerably depending the lines studied 
and the experimental conditions. obtain the strongest components the 


TABLE 


Previous work | Previous work 
Present work AA 1} Present work XA 
Transition A Tre\* vyae, Werner Schiiler Transition r, A Tre vvace Cm~! Werner Schiiler 
3)P-4'1D 4788.363 20878.13 8.8 8.5 || 33D-6F 2728.288 5 36642.17 8.4 
2934.022 


Nore: tDiffuse line; sh. ‘‘shoulder’’, incompletely resolved weaker companion of a stronger line. 
“The estimated intensities in this column are relative intensities within each multiplet. The numbers given for 
different multiplets are not comparable. 

‘Series and Willis (1958). 

‘Note added in proof: In his thesis, which has only now become available to us, Werner (1927) gives the 
following additional lines: 2 'S-3'P 1420.7 A, 2'P-#!'S 1253.8 A, 2 4S 1166.4 A, and 2 1D 
1132.1 A. 


Wavelengths and wave numbers (except 5485 see Table 
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Sections the spectrum the vacuum ultraviolet. 
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Fic. Structure the line 5485 (a) for natural lithium, and for lithium 
which was enriched each case two exposures different density are given. The 
components are numbered Table 
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line 5485 the second order the 21-ft grating only seconds were 
required, while some the higher members the series required 
hours. 

The plates were measured means photoelectric comparator similar 
that described Tomkins and Fred (1951) and built Stoicheff 
this laboratory. Most lines were measured several plates. The average 
wavelengths and wave numbers all observed lines except the line 
5485 are given Table and compared with the previous measurements 
Werner (1925, 1926) and Schiiler (1924, 1925, 1926). The infrared line 
was only recently observed under low dispersion Series and Willis (1958). 
comparison measurements different plates estimated that the 
measurements above 2000 are accurate +0.003 while those below 
2000 are accurate +0.001 This corresponds throughout wave 
number accuracy better. Estimated relative intensities are 
given only for triplets and apply only within each triplet. 

The line well known have somewhat complex structure 
(Schiiler 1925, 1927, Granath 1932). Figure shows this line obtained 
the second order our 21-ft grating. This spectrum was obtained with 
natural lithium Li’, were most the other spectra. Figure 
shows spectrum obtained with sample which the isotope 
was enriched 96%. Comparison these two spectra makes possible 
sort out the Li® and Li’ lines. These are listed separately Table and 
compared with the data Schiiler (1927). may noted that our wave- 
lengths for the components this line are all about 0.21 than 


TABLE 


Fine structure and hyperfine structure the line 5485 Li? and 


air, obs. Obs. 
Component (present Triplet Hypertine 
number Present work Schiller Jest work) comp. comp. Yvaec. Cale Teal 
Li’ (1) 54838 3458 3.55 1 18231 .975 1/2-3/2 a 0.520 
(2) 5483 4590 3.67 2 18231.598 } 3/2-3/2 1.179 
(3) 5483 3. 18230.94, 5/2-3/2 2.337 
1/2-3/2 d 1.191 
1/2-1/2 f 1.541 
§/2-7/2 g 8.000 
5 5 85 ) 8228.18: 28S:-23 
(5) ABS: 4.6) 1 1 P: \ 3/2-3/2 h 448 
(6) 5484. 5059 1.785 5 18227 .8l» 5/2-8/2 18227.8 2 713 
O97» 1 18227. 480 $/2-3/2 18227. 48: 0.361 
1/2-3/2 18226.57,4 2.289 
(8) U7 36 5.17% 18226. 56, | 4 3/2-3/2 m 18226. 525 2.713 
3/2-3/2 o 18226. 18s 0.37 
(13) 5485. 5.31 5 18226 2°S,-2 * Pi 32-12 18226050 
5485. 1767 5.376 7 18225. 5/2-3/2 is 3. 287 
(10) 5485. 2805 5.47 3 8225. $/2-3/2 1.337 
Lit (14a) 5483. 4.10 3 18230. 199 s 18230. 19 
(14d) SASS 18230. Otic 2-1 18230. 06, 


(il) 5ABS 5.05 5 18224. O85 v 95s 
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Schiiler’s. This shift due the fact that Schiiler measured second order 
against the line 2741.2 appearing fourth order assuming old value 
2741.29 for this line which more recently (France 1930; Datta and Bose 1935; 
Johansson 1959) has been revised 2741.184 This change accounts for 
the difference. Our wavelengths, based measurements five different 
plates, are believed supersede those Schiiler. 


TERM SYSTEM, FINE STRUCTURE, AND HYPERFINE STRUCTURE 

The assignment most the observed lines has already been made 
Werner (1925, 1926) and Schiiler (1924, 1926). Fig. energy level 
diagram plotted showing the observed lines. The only unassigned lines 
multiplets are 3714.2 and 2934.2 already observed Schiiler (1925) and 
Werner (1926) and 4607.4 and 2337.0A previously observed Werner 
(1926). These lines not result from transitions between levels shown 
the diagram Fig. difficult certain that they are due Lit although 
they behave the same way the Lit lines. they are due Lit they must 
ascribed transitions between doubly excited states. 


(20,000 | 2 


140,000 } 


Fic. Energy diagram observed levels except the L4S ground state are 
shown and the observed transitions are indicated. (See, however, footnote Table 
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With the exception the lines and all the classified 
lines observed the visible and near ultraviolet regions have their lower 
states three-quantum states. All transitions between the three-quantum and 
two-quantum states lie the vacuum region. 

The detailed analysis the spectrum and particularly the determination 
the series limits made complicated the fact that the hyperfine structure 
splittings are the same order even larger than the ordinary fine structure 
splittings. This complication has been well known for many years for the line 
5485 corresponding the transition But this the only line 
whose fine structure had previously been resolved (Schiiler 1925, 1927, 1930). 
well known the nuclear spins Li? and are and respectively. The 
present analysis based these values and way may said confirm 
them that were necessary. 

shown Table II, slightly better resolution has been obtained for 
the line than Schiiler. His assignment the Li® lines has been 
confirmed the observation the spectrum nearly pure Li® (Fig. 
this spectrum, component (14) has been partially resolved into two components 
which must correspond the two stronger the three expected components 
(0-1, 1-1, 2-1). Indeed the separation these two components agrees well 
with that predicted from the known ratio the magnetic moments 
and Li® 3.9614, see Kopfermann (1958)] and the observed 
splitting factor 0.257) the component Li’. This pre- 
dicted value 0.130 while observed. The lines (7) and (11) Li® 
are not further resolved since the hyperfine splitting and 
the same order that and there greater number line components 
within the same wave number interval. 

Taking account the statistical weights, the center the state 
Li® the ordinary mean the levels and and therefore the 
mean the components and when subtracted from components 
(7) and (11) gives the positions and below One finds 3.07, 
and i.e., well known, the state partially inverted (see 
Bethe and Salpeter 1957). With the aid these triplet intervals can 
now obtain calculated set hyperfine levels both for Li? and For Li® 
the standard hyperfine structure formulae (see Kopfermann 1958) 
used since the hyperfine splitting small compared the triplet splitting. 
However for Li’, since the hyperfine splitting much larger, the more elaborate 
method Giittinger and Pauli (1931) has used which takes account 
the perturbations between levels the same 

Giittinger and Pauli have used the same splitting factor for both 
and While indeed the main contribution the splitting made the 
expected make small contribution (calculated 0.019 from formula 
(26.9) Kopfermann) while the contribution the electron should 
entirely negligible. Therefore the factor for expected somewhat 
smaller than for 24S. Accepting the above predicted contribution 0.019 
the electron, one would expect 0.495 for 0.514 for (as 
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obtained from the three components 3). However, find that 0.495 
for gives slightly poorer representation the observed components 
than 0.514. Thus adopt (in agreement with Pauli and Giittinger) the 
same splitting factor for 2*P and The energy levels 
derived from this splitting factor and the observed line components are listed 
Table III for both Li? and The wave numbers the lines derived 


TABLE III 
Energy levels and Li? and Li® giving the best fit with the 
observed fine structure components the line 5485 


| 
| 
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Li? 
Multiplet energy above center energy above center 
7/2 18228 604 


from the energy levels Table III are included the second last column 
Table The agreement with the observed wave numbers, particularly 
unblended components, satisfactory. The revised triplet splittings are: 


The calculated intensities are given the last column Table These 
were obtained the method Giittinger and Pauli but without some 
the approximations made them. The agreement with experiment slightly 
better than found Giittinger and Pauli; for example, the intensity ratio 
components (3) and (2) was observed Schiiler (1930) 1.95+0.04. 
According Giittinger and Pauli this ratio was 1.51 and this dis- 
crepancy was considered serious Schiiler. However, the revised 
calculations give ratio 1.98, excellent agreement with experiment. 

None the other lines Lit has been resolved completely but some 
resolution has been obtained for all triplet transitions and few cases 
even for singlets. there were hyperfine structure one would expect only 
and should single because the triplet splitting 
for with and for all and states would expected 
below 0.2 i.e. unobservable with the resolving power our gratings. 
which could formally explained due splitting but five 
times large the expected triplet splitting. addition, each the singlet 
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lines has satellite appreciable intensity even though the basis 
ordinary hyperfine structure theory splitting singlet levels expected. 

order understand these apparent anomalies necessary remember 
that the hyperfine splitting and similar two-electron systems deter- 
mined mainly the electron, present all levels considered, and that 
therefore the splitting all triplet levels has about the same magnitude 
the state (see Table Even the singlet levels will show such 
splitting their separation from the corresponding triplet levels (of equal 
smaller than the hyperfine splitting the latter. Unfortunately not 
possible use standard hyperfine structure formulae since the coupling 
conditions are quite different from the usual ones. While for small hyperfine 
splitting the coupling may written symbolically 


(S82) = (LS )I = JT = F 
for large hyperfine splitting and small triplet splitting have (with 


may even happen that the hyperfine splitting large compared with the 
energy difference. this case the coupling may described 

case loses its meaning. Instead there quantum number corre- 
sponding the resultant and which very loosely coupled. 
case the nuclear spin couples with (or ji, form resultant 
F,, the value the ion; very slightly coupled with and before 
not defined. 

Case can treated the method which takes into 
account the interaction levels the same applying this method 
state with small triplet splitting, one finds indeed that the resulting levels 
are arranged three groups corresponding toG 1/2, 3/2, 5/2 and for each 
for 3°D and Fig. (Such arrangement can also seen the 
levels discussed Fred, Tomkins, Brody, and Hamermesh (1951)). 
The splitting between levels the same decreases rapidly with increasing 
and 

Since pure spin quantum number (electron spin nuclear spin) one 
would expect the selection rule hold. result the transition 
between two states belonging case should consist one strong line 
(AG 0), possibly with weak satellites (with 0). The observed structure 
the lines (Table might have been understood this basis: 
there strong central line (AG 0), very weak short-wave satellite 
which would fit the transition from 5/2 toG 3/2 (i.e. 
+1), and much stronger and closer long-wave satellite 0.25 which 
would correspond the splitting the sublevels with 3/2 the 
lower state, e.g. the transition from 7/2 5/2 (see Fig. with 
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(a) (b) 

(b) using case (1) for and case for and disregarding all smaller interactions. 

However, this interpretation cannot maintained because additional 
complication arises from the fact that the levels have very nearly the 
same energy 'F, that is, that the energy difference much smaller 
than the hyperfine structure splitting. Therefore coupling scheme (II) must 
used for the states. first approximation the splitting determined 
that the ground state the ion which 0.995 (calculated 
from the known magnetic moments and and the observed splitting 
hydrogen means standard formulae, see Kopfermann (1958)). Thus 
for large expect two component levels with and separated 
approximately With decreasing each these levels would 
expected split into several component levels for and for 
even for the lowest levels Since not defined the two main 
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would expect both them combine with the state. now very 
significant that the stronger members the series are indeed observed 
double consisting one strong and one weak component which are 
separated 1.07, 1.00, and for and respectively (see 
Fig. 3). Within the accuracy the measurements this splitting approaches 
the predicted value 0.995 the ground state. must noted 
that the stronger component the shortward side, that is, corresponds 
the upper component level with 

very similar anomaly has been observed previously for the 
series Paschen (1923) and Sawyer and Paschen (1927), and has 
been interpreted Paschen (1932), Goudsmit and Bacher (1933), Paschen 
and Ritschl (1933), and more recently and greater detail Suwa (1950, 
1951). The latter authors have also discussed the effect this peculiar coupling 
the corresponding triplet transition Similar transitions 
have been observed and discussed Paschen (1935) and Paschen and 
Campbell (1938). 

While Al* and the lower state shows fairly normal triplet 
structure, for the lower state the triplet transitions has very small 
triplet splitting which completely overshadowed the hyperfine splitting 
leading coupling scheme (1). Thus have transition between state 
*F) belonging coupling scheme and state belonging 
coupling scheme (1). Rigorous intensity rules for such case have not been 
worked out but may attempt apply the same correspondence con- 
siderations that were used Goudsmit and Bacher, namely that those 
transitions which the relative orientation the spin the electron 
and the nuclear spin remains unchanged are the strongest. Fig. 
energy level diagram for the transition shown neglecting all the 
weaker interactions, i.e. there are only two levels the upper state (F; 
and three the lower state 1/2, 3/2, 5/2). addition the level 
also shown. The range values indicated for each level. The selection 
rule for allows six transitions, for and two, and 
lower state the spins and are parallel, and 1/2 they are 
antiparallel (see the arrows‘to the right the levels). Therefore according 
intermediate, and and forbidden. can seen comparing Fig. 
with the spectrum, Fig. and Table the strongest one the three observed 
components can interpreted the superposition the two transitions 
and while the other components can interpreted and 

close enough each other for coupling (II) considered except very 
large Therefore all transitions other than can dealt 
with the basis coupling scheme that is, application the 
method. 

The fine structure and hyperfine structure the transition can 
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means elementary fine structure theory (see Bethe and Salpeter (1957)) 
and Giittinger and Pauli’s method. The triplet splitting for predicted 
figures are and respectively). this basis, relative 
positions the component levels can derived. Fig. the calcu- 
lated structures and are compared. Corre- 
sponding transitions are designated the same letters (see Table 
seen that accordance with the selection rule which holds coupling 


83467 83465 


27133 27130 


Fic. Theoretical splitting patterns (a) (b) (c) Note 
that the wave number scale (a) half that (b) and The lines marked correspond 


case with increasing the intensity more and more concentrated the 
rapidly lose their intensity. the same time the transitions with 
move closer and closer together. With sufficient resolution the line 
would consist three components with spacings 0.61 and 0.32 However, 
the vacuum ultraviolet such resolution difficult obtain and the 
present work only one line has been observed. From the position this line 
the energies the levels can derived given Table IV. 
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For the structure represented Fig. 8(c) expected. Because 
the smaller triplet splitting, the lines with have still lower intensity 
than for and the strong lines with are more concentrated 
near the center. agreement with this prediction only single line has 
been found. Even with the higher resolution the near ultraviolet the close 
components near the center cannot resolved. Similar considerations apply 

The transition expected have structure very similar that 
except for the slight change hyperfine splitting compared 
Since this line the vacuum region complete resolution cannot 
expected, but three components are indeed found shown Fig. 
Their spacings are observed 2.02 and These spacings agree 
well with the spacings the centers gravity the corresponding predicted 
groups lines (see Fig. which are 2.03 and 3.08 

Again the near ultraviolet series similar those just discussed. 
fact the structure the line except for the very slight change the 
hyperfine splitting (and the reversal the scale), should identical 
with that shown Fig. Three components 
observed whose spacings (0.44 and 0.27 agree reasonably well with the 
spacings the centers the three expected central line groups (0.63 and 0.34 
The spacing appreciably different from the predicted triplet spacings 
(0.904 and used the computation Fig. The strong 
effect the perturbations the hyperfine structure (affecting both position 
and intensity the lines) clearly shown this way. 

For the transition rather complex structure expected since 
the state still far removed from coupling case (I) while approach- 
ing it. The selection rule longer valid and instead the 
transitions possible when the hyperfine structure small compared the 
fine structure, here all transitions compatible with the rule 
can occur. Figure shows the intensity distribution derived from detailed 
calculation according the method. Unfortunately the 
available resolution the vacuum region allows the observation only three 
components (see Fig. whose spacings (2.85 and 1.88 agree satis- 
factorily with those predicted from Fig. averaging the three line 
groups (2.94 and 

The transitions the visible region are similar but 
the lower state closer coupling case (1) and the upper state almost 
perfect case (1). Therefore the structure somewhat simpler that the 
transitions with are more strongly emphasized compared the other 
transitions. Figure shows the calculated structure for Fig. 
the main transitions are shown energy level diagram. the lower 
state were good approximation case (I) the upper state, only 
one line would observed which would the superposition the three 
group (a) three lines the left Fig. 10. Because the deviation from 
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Fic. Energy level diagram for showing principal transitions. 
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case the lower state splitting arises even for the transitions with 
the line (c) arises from The three observed lines correspond apparently 
the three groups and (c) with few transitions mixed 
(see Fig. The observed spacings 0.28 and 0.39 agree satisfactorily with 
those predicted from Fig. i.e. 0.303 and and the estimated 
intensities 5:5:1 are compatible with the predicted 15:11:5. 

the basis the preceding discussion possible establish the relative 
positions the various energy levels the triplet system. For the singlet 
system the situation very much simpler and need not 
detail. The relative position the singlet and triplet system determined 
the singlet and triplet Bergman series since for both the upper states are 
the same, i.e. and (see Fig. 7). Tabel the energy values all 
observed levels above the ground state are listed including, where possible, 
the individual hyperfine levels. the case the and and levels, 
the calculated energy values are listed that give the best representation 
the observed incompletely resolved hyperfine structures. The absolute accuracy 
the energy values since they are based approximate 
values the far ultraviolet lines (see below) but the relative accuracy the 
energy values for most excited states should except that the 
levels have greater possible error relative 
the other levels (perhaps +0.15 since they are related them only 
series limits. 


SERIES LIMITS, TERM VALUES, AND IONIZATION POTENTIALS 


the usual methods determining series limits (plotting against 
are applied the observed three singlet series and four triplet series (using 
the strongest lines the latter) the limits listed Table are obtained. The 


TABLE 
Series limits 


designation Limit, values used 


1 


“Weak satellite. 
Strongest component. 
“Average of two strongest components. 


two series and 'D, expected, give very nearly the same limit, 
48327.20 and 48326.95 respectively. The latter value appears somewhat 
more reliable since the magnitude the Rydberg correction (and therefore 


Series 
_ ‘ 
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also the higher terms the Ritz formula) much smaller for than 
for must, however, add correction account the splitting 
the ground state Since for all observed values the and states 
are above the corresponding triplet states they must correlate with the 
upper (F; component the ground state. Referring the center 
gravity the ground state, i.e. the energy that would have 
without the effect nuclear spin, have subtract 3/8 the splitting, 
that 0.37 and therefore obtain for the term value 


48326.58+0.10 


Similarly the limit the Bergman series (strong component) corresponds 
the upper component the ground state discussed detail the 
preceding section. Correcting the unsplit ground state obtain 


and combining with the wave numbers the lines and 
108270.81 


deriving term values from the limits the triplet series, the fine structure 
and hyperfine structure must taken into account. According the inter- 
pretation Section the observed limit the strong component the 
Bergman series gives the term value the 5/2 component 
Fig. Referring again the unsplit ground state the ion 
find 


48835.07 
48836.08 


The strongest observed component the transition corresponds 


The term values the other components are readily obtained from 
this number and the energy values Table III. 

From the same table obtain for the lowest triplet state (weighted average 
value can obtained from the series. Only the first two members 
this series are the two components and Fig. resolved. The series 
limit corresponds the average these. The 5/2 group levels 
n*D must the level the ion and therefore according Fig. 
the series limit gives the average the levels and below 
below Adding the average energy and above 
the one hand and that the other 


however, added the end this paper. 
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(Table IV) and correcting the unsplit ground state the ion yields the 
values 134044.24 and The mean these two 
values, agrees satisfactorily with the value from the Bergman 
series but perhaps somewhat less accurate because the more complicated 
hyperfine structure and the greater Rydberg correction for than for 
states. may consider the mean the final 
value based triplets only. However, since the relative positions the 
singlet and triplet levels fixed (see Section can obtain also 
from given above using the separation from Table and 
find 134044.16 The agreement the two values very 
satisfactory and gives confidence the determination the series limits. The 
last value probably somewhat more accurate since based single 
lines. give twice the weight forming the final average which 


This must considered present the best available value. may 
compared with the best previous value, 134033 given Moore (1949). 
While have not yet remeasured the far ultraviolet resonance lines 
may use the wave numbers given Edlén quoted Robinson (1937) 
for together with our new term values for and 
4'P order obtain revised value for the ionization potential Lit. 


find 


Using weights corresponding Robinson’s quoted error limits obtain 
final mean 


appears that the major part the previously quoted error (+25 
was due the uncertainty the exact values the term values the 
excited states. The error quoted Robinson for the wave number the 
resonance lines was and this presumably the limit error the 
ionization potential. The excellent agreement the three independent values 
suggests that the probable error the mean may less, perhaps 
Indeed private communication received after the preceding discussion 
was written, Edlén has pointed out that the error estimate the wave- 
lengths the resonance lines was conservative one. also emphasized 
that selecting standards for these measurements use was made 
hydrogen-like spectra calculated without knowledge the Lamb shift, 
error that does affect several other wavelength determinations the region 
below 200 
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CONCLUSION 


Hart and Herzberg (1957) calculated the energy the ground state Lit 
means 20-parameter Hylleraas function and obtained (including 


610085.6 


These calculations were much improved Pekeris (1958), who went 
203 parameters with certain approximations and obtained 


610087.44 


value that estimated accurate Subtracting this 
Dirac value from the observed value, one obtains experimental value 
for the Lamb shift 


with error possibly The predicted Lamb shift —7.8 
(Pekeris 1958). The agreement better than might have been expected 
the basis the error the observed value. 

For the state until now only 27-parameter non-relativistic energy value 
has been calculated (Hart and Herzberg 1959) 


The difference this value and the observed value, the sum 
the relativistic correction and the Lamb shift. Once the relativistic correction 
has been calculated observed Lamb shift with accuracy 
will obtainable. 

While the present paper the excited energy levels have been established 
with adequate accuracy, definitive determination the ground state 
remains done. For this purpose the far ultraviolet resonance lines will 
have measured with accuracy better than 0.001 far the 
present results go, they confirm the predictions quantum electrodynamics, 


did those He. 


NOTE ADDED PROOF 

Professor Edlén has called our attention the fact that the Rydberg 
correction the state does not fit very well with the Rydberg corrections 
the higher states and that Werner (1927) his full paper (not accessible 
the time writing our paper) has given different value for the 
line was not mentioned Werner’s Notes Nature and was not observed 
our work. Werner’s identification this line correct, would imply 
that the infrared line found independently Series and Willis 
(1958) and ourselves incorrectly assigned and that the correct 
line must 9572 The search for the infrared line, both Series and 
Willis’ work and our own, was based extrapolated value for the 
state given Werner (1926) and quoted Moore (1949). did not search 
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near 9572 because the much lower sensitivity the infrared plates 
that region. While the Rydberg correction the state based our 
identification indeed somewhat large, the Rydberg correction based 
identification the 1420 line somewhat small. the present 
moment difficult decide which value the correct one. even 
possible that neither one the two values correct. None the remaining 
conclusions our paper affected this uncertainty the state. 
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PRINCIPE ANALYSE DES COURBES EXPERIMENTALES 
COMPOSEES: APPLICATION MESURE DENSITE 
DES GRAINS DANS LES TRACES FORTEMENT 
IONISEES DANS LES EMULSIONS NUCLEAIRES! 


RESUME 


expose principe d’une méthode mathématique, utilisant les transformées 
Fourier, qui rend possible décomposition, dans ses éléments constitutifs, 
courbe expérimentale complexe, lorsque celle-ci provient superposition 
plusieurs courbes ayant toutes méme forme analytique, mais des amplitudes 
considére, particulier, possibilité d’utiliser cette analyse 
pour déterminer, partir des courbes photométriques, densité des grains des 
traces laissées par les particules fortement ionisantes dans les émulsions iono- 
graphiques. montre dans exemple simple, idéalisé, 
précision méthode. 


INTRODUCTION 


Dans présent compte rendu expose principe d’une méthode d’analyse 
des courbes expérimentales complexes lorsque celles-ci présentent des carac- 
mathématiques simples. Cette méthode peut étre utilisée, général, 
dans tout probléme physique, lorsqu’il s’agit décomposer dans ses éléments 
constitutifs, une courbe qui provient superposition plusieurs courbes 
ayant toutes méme forme, mais des hauteurs différentes. 

premier exemple d’application cette méthode est celui 
des courbes photométriques représentant des trajectoires des parti- 
cules ionisantes dans une émulsion ionographique. 

Cette analyse pour but déterminer densité des grains développés 
leur opacité degré noircissement dans partie finale trajectoire 
lorsque les mesures visuelles deviennent imprécises méme impossibles. 

mode d’analyse été pour s’appliquer des 
grains, dans d’une ‘‘famille simplement 
(Ahmad Demers 1959). Dans cette hypothése, les courbes photométriques 
des grains ont méme forme analytique une constante multiplicative prés; 
plus, lorsque les grains sont suffisamment rapprochés, les courbes photo- 
métriques montrent une région superposition sans interférence, avec simple 
addition des densités optiques. Nous renvoyons ces auteurs 
pour justification ces hypothéses, qui paraissent rendre compte fagon 
trés convenable des observations photométriques: absence d’interférences 
perceptibles possibilité d’appliquer modéle méme pour des grains écartés 
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FORMULATION METHODE 
s’agit d’analyser une courbe qui est mathématiquement exprimée par 


relation: 


Dans cas physique actuellement sous considération, F(x) représente 
trace selon son axe sont des courbes méme famille dont 
l’allure représente des grains individuels long tra- 
jectoire. Ces courbes sont symétriques par rapport aux points qui sont 
les centres des grains leur forme analytique est arbitraire, mais 
supposée connue dans les applications; une bonne approximation semble 


étre 


avec Définissons une fonction distribution des centres des grains 


Alors, cause des propriétés des fonctions Dirac, 


limite supérieure est prise par excés sans introduire par cela 


aucune erreur. 
prenant transformée cosinus Fourier chaque membre 


(4) tenant compte théoréme des convolutions des fonctions trans- 
formées Fourier (Morse Feshbach 1953), obtient relation simple: 
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(8) D(t) g(t) cos dt. 


pourra donc obtenir distribution des grains D(é) partir trans- 
formée cosinus quotient des transformées cosinus des courbes photo- 
métriques d’une portion trace F(x) d’un grain individuel, 

est, dans les hypothéses admises départ, une somme fonctions 
représentables physiquement par une série traits verticaux placés sur 
chacun représente centre d’un grain, hauteur ces traits 
prise égale coefficient signifiant l’opacité des grains question. 

suite mathématique, fonction expérimentale continue 
F(x) trouve remplacée par une fonction discontinue 
approximation, qui est probable, fonction 
résultante sera certainement modifiée. est difficile d’imaginer d’avance 
ces modifications sans calculs numériques précis; peut espérer que 
courbe photométrique trouvera alors remplacée, sinon par une série 
traits infiniment étroits, moins par une nouvelle courbe continue 
position des grains seront plus facile reconnaitre. 


DISCUSSION 


Dans présent compte rendu, exposera seulement quelques considéra- 
tions sur résolution mathématique probléme. 

d’abord déterminer transformée qui sera méme pour 
toutes les traces des plaques d’une méme émulsion; attribuant une 
forme analytique convenable, pourra facilement calculer obtenir des 
tables, s(t). 

calcul f(t) est difficile car F(x) peut pas étre approchée par une 
fonction analytique. Généralement faudra utiliser calculateur électronique; 
dans certains cas, cependant, pourrait songer développement 
Legendre: 


Une fois obtenus les coefficients calcule aisément transformée 


avec les fonctions sphériques Bessel habituelles. 

est difficile prévoir g(t) sans effectuer les calculs; peut 
affirmer cependant devrait pas étre singuliére car, général, 

Par conséquent, seule difficulté théorique qui apparait provient fait 
que est une fonction Dirac qui rend problématique 
calculateur électronique. Cependant calcul numérique est encore 
possible car, grace limité, obtient une fonction 


PASKIEVICI: COURBES EXPERIMENTALES COMPOSEES 1317 


approchée celle Dirac qui tend vers celle-ci lorsque 


augmente indéfiniment. 


possible connaitre précision des résultats obtenus, commencé 
programmation d’une série calculs partant des valeurs numériques 
d’une courbe construite synthétiquement propos. 


APPLICATION NUMERIQUE 
courbe synthétique (figure été composée huit courbes individuelles 


ayant forme 


hors cet intervalle. 


coefficient dépend longueur utilisée; dans notre cas, 
donne les valeurs dans méme unité, unités arbitraires; 
choix ces valeurs correspond point vue physique une trace 
ionisante complétement saturée, lorsque comptage des grains devient 


impossible. 


TABLEAU 


courbe F(x) été calculée par points, pour valeurs, utilisant des 
tables quatres décimales exactes; voit qu’elle présente uniquement deux 
maxima. Les calculs numériques ont été effectués d’un calculateur 
arithmétique Librascope LGP-30 Royal McBee Corporation; pour évaluer 
les intégrales, utilisé méthode Simpson. 


20 60 ) 


1318 CANADIAN JOURNAL PHYSICS. VOL. 37, 1959 


Les figures représentent, respectivement, les transformées 
s(t) leur quotient g(t); dans derniére figure tracé aussi, pointillé, 


f(r) 


g(t) 


a 


fonction g(t) calculée exactement partir fonction connue dans notre 
cas, dont elle est transformée. Plus explicitement 


comparant ces deux derniéres courbes, constate une excellente pré- 
cision jusque vers 40, qui correspond premiére moitié 
valle considérée, tandis qu’elle diminue légérement dans seconde moitié, 
comme fallait attendre d’ailleurs, cause croissante sur 


O(x) 
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g(t) définie jusqu’a égal respectivement. Dans cette derniére 
figure, dans laquelle indiqué aussi, pointillé, transformée Fourier 
g(t) exacte, observe aisément les huit pics principaux qui correspondent 
aux huit courbes initiales. L’analyse courbe F(x) est ainsi réalisée. 

est important remarquer que courbe figure devrait étre 
mais qu’elle n’est qu’une approximation d’une somme fonctions Dirac 
ceci est fait que d’intégration sur est limité une 
certaine valeur obtient pour 


sin (E—E£;)to 

est suffisamment grand, comme dans notre cas 80, séparation 
entre les courbes est suffisante pour déterminer avec précision nombre 
maxima leur emplacement; les valeurs des maxima qui sont proportionnelles 
aux coefficients d’opacité sont pas déterminées avec méme précision 
parce que les courbes individuelles présentent encore des domaines d’inter- 
férence. 


CONCLUSION 

méthode développée ci-dessus devrait permettre, dans 
décomposer une courbe, qui est une superposition plusieurs courbes 
individuelles, dans ces éléments constitutifs. Les résultats obtenus partir 
d’une courbe synthétique, mais rapprochant autant que possible 
courbe réelle, montrent que décomposition est mathématiquement possible. 

Une application immédiate est vue: celle d’obtenir position par 
suite densité des grains leur opacité, long des traces ionographiques 
forte densité, lorsque comptage habituel devient impossible. 

des profils expérimentaux présente des 
dues essentiellement caractére approché base qui con- 
sidére tous les grains méme parfaitement alignés selon axe, 
aux erreurs expérimentales. L’influence ces facteurs dans procédé 
présentement développé sera étudiée dans une publication ultérieure. 

connaissance ces valeurs contribuera structure des 
traces, perte par ionisation, discrimination plus rapide 
plus précise des particules; elle pourra vérifier, particulier, 
certaines hypothéses avancées sur des grains 
contagion dans ionographique. 
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dans ses derniéres mesures. 
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NOTES 


SOME OBSERVATIONS NUCLEAR ALIGNMENT COBALT METAL 


The possibility obtaining nuclear alignment using the interactions 
ferromagnetic single crystal was demonstrated Grace al. (1955, 1957), 
who showed that when cobalt metal single crystal containing little 
cooled about 0.05° spontaneous nuclear alignment occurs parallel 
the hexagonal axis the cobalt crystal. This alignment causes the y-rays 
from the emitted anisotropically with respect the hexagonal 
axis the cobalt crystal. The mechanism which produces this alignment can 
described, qualitatively, follows. Cobalt metal crystallizes the hexa- 
gonal system; also ferromagnetic and the direction spontaneous mag- 
netization the domains parallel the hexagonal axis below 250° Thus, 
low temperatures, all the unpaired electron spins responsible for the ferro- 
magnetism are aligned parallel the hexagonal axis, and each spin produces 
magnetic field several hundred thousand gauss, parallel this axis, 
the site the nucleus which attached. temperatures below 
the cobalt nuclei are aligned parallel the hexagonal axis, because this 
orientation describes the state lowest energy, which preferentially popu- 
lated. have aligned cobalt metal single crystal this method, 
and have made measurements the anisotropy y-rays emitted the 
aligned nuclei. Our results differ significantly from those other 
workers, and account them presented below, with possible explana- 


tion. 

piece cobalt metal single crystal, size, was 
given Dr. Myers, formerly the Department Mining and 
Metallurgy, University British Columbia. was irradiated activity 
about the neutron flux N.R.X., Chalk River. was cooled 
temperatures the order 0.04° being placed contact with pieces 
potassium chromium alum, cooled this temperature adiabatic demag- 
netization. The temperature this potassium chromium alum was measured 
the ballistic method. The intensity emission y-rays was measured 
scintillation counters placed along, and perpendicular to, the hexagonal axis 
the cobalt crystal. This apparatus has been described greater detail 
Daniels (1958). Various methods were used ensure good heat 
contact between the cobalt metal and the chromium potassium alum. These 
were: 

(i) The crystal was placed between two cylinders, in. diameter in. 
long, compressed chrome alum, and the gap between was filled with glyc- 
erine, with solution chrome alum 50% water, with sludge 
Apiezon oil and chrome 
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(ii) The crystal was pushed into hole cut copper strip 3cmX1 
and the whole assembly was copper-plated, bonding the crystal the copper 
fin. This was sandwiched between four half cylinders chrome alum, and 
lubricated with sludge chrome alum and glycerine. 

(iii) About 200 No. B&S copper wires providing surface area 
about 150 cm? were soldered the fin. Chrome alum powder moistened with 
Apiezon oil was carefully placed between the wires, and the whole assembly 
compressed into solid pill press. 

(iv) thin copper foil about 200 cm? surface was soldered the copper 
fin, and folded concertina shape. Chrome alum powder and Apiezon oil 
was placed between the folds and pushed down hand-tight. The assembly 
was held lucite container. 

The results obtained from all these methods were self-consistent during the 
time that the anisotropies were decreasing and the chrome alum was warming. 
Hence are satisfied that there was adequate heat contact between the 
chrome alum and the cobalt crystal all our experiments. 

The anisotropy emission y-rays, defined 


where and are respectively the counting rates along and 
pendicular the hexagonal axis the cobalt crystal. The results typical 
experiment are given the figure. assume that after about hour 
the cobalt crystal has come into thermal equilibrium with the chrome alum, 


0.045°K 
€ 12 
0.08 
0.04 


Fic. The measured anisotropy y-ray emission from cobalt cobalt metal single 
crystal, and the temperature the specimen, both function time. This for our 
cobalt crystal, after heat treatment, cooled surface No. B&S copper wires 
described the text. The solid points represent the values obtained Oxford for the 
temperatures indicated. 
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can correlate with the absolute temperature high temperatures 
(above about 0.05° this case) where constant which 
describes, among other things, the hyperfine coupling the nucleus 
with the electron spins. 

observed This considerably larger than the Oxford 
values, (Kurti 1959), and rather smaller than the value 
reported Khutsishvili (1955). Now the hyperfine coupling 
responsible not only for the anisotropy y-rays, but also for the specific 
heat below This has been measured Heer and Erickson (1957), Heer 
(1958), and Arp al. (1957). Heer and Erickson give for the specific heat 
hyperfine coupling the form have 

= (A /12 and 39 (A 60 — Boo) / 28k°7 
where the subscripts and refer the isotopes and re- 
spectively. The term negligible, hence can calculate what value 
should expect from the measured values specific heat. The two values 
quoted give, respectively, and 

The crystal was heated red heat for few minutes air and the 
experiments repeated. This time the anisotropies were much smaller, 1.5 
suggest the following explanation for the discrepancy between all 
the values quoted for Cobalt has phase transformation 417° the 
low-temperature stable phase being hexagonal, and the high-temperature 
phase cubic. result, almost impossible prepare single crystals 
hexagonal cobalt, small cubic regions being present stacking faults. 
proportion 10% cubic phase quite common (Houska and Averbach 
1958), and proportions more than 50% have been observed supposedly 
hexagonal cobalt (Edwards and Lipson 1942). assume, seems reason- 
able, that our crystal contained more the hexagonal phase before heat 
treatment than after, then must assume that the hyperfine coupling 
the cubic phase much less than the hexagonal phase. This would explain 
why different experimenters have obtained different results for and also 
why the different measurements specific heat cobalt metal differ 
amounts which are apparently much greater than the experimental error. 
would appear that the crystal used Khutsishvili contained more the 
hexagonal phase than did ours, and that ours contained more the hexagonal 
phase than did the crystal used Oxford. would interest measure 
both the specific heat and the anisotropy y-rays from cobalt metal for 
which the proportions cubic and hexagonal phases are known. The pro- 
portions can measured X-ray analysis, however, have not the 
necessary equipment for this our disposal. 


The support for this work the National Research Council and Atomic 
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NUCLEAR SPIN RELAXATION LOW-TEMPERATURE LIQUIDS 


The nuclear spin lattice relaxation time may expressed terms 
the correlation functions those interactions which enable energy 
exchanged between the nuclear spin system and the other degrees freedom 
the material which the spins are located (Bloembergen, Purcell, and 
Pound 1948). Measurements are therefore capable giving information 
related the statistical mechanical properties matter. this note 
discuss the interpretation measurements liquids, and particularly 
low-temperature liquids, where the potentialities the nuclear magnetic 
resonance technique have not yet been realized. 

There are two main types relaxation mechanisms considered. 

the cases atomic liquids, such helium, xenon, etc., the spin 
relaxation caused the dipole-dipole interaction between neighboring 
nuclear spins and the well-known result for expressed terms the 
correlation functions where are the polar co-ordinates 
the vector joining pair atoms (Bloembergen, Purcell, and Pound 
1948). 

the molecules have rotational angular momentum additional mechanism 
present due interactions within the molecules which couple the nuclear 
spins Here the problem calculate the correlation function 
(J(t)J(t+7)). The simplest case low temperatures where all the 
molecules are the state and may written 
terms the mean time spent molecule state, which turn 
depends the correlation functions the anisotropic part the inter- 
molecular potential (Bloembergen, Purcell, and Pound 1948; Moriya and 
Motizuki 1957; Moriya 1957). 


*National Bureau Standards, Washington, D.C. 
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will shown future publication, for many systems plausible 
approximation for the correlation function which 
function only the vector joining the pairs and not their momenta, 


where P(R, R’, related the probability per unit volume that two 
Note that the form which the radial distribution function g(R) appears 
equation (1) slightly different from the usual form which appears 
the correlation function (Torrey 1958). For liquid assume the particles 
obey the diffusion equation (Chandrasekhar 1943) which case 


where the diffusion constant. 

g(R) has been evaluated using the superposition approximation for liquids 
which molecules interact via Lennard—Jones potential (Kirkwood, 
Lewinson, and Alder 1952). applying these results low-temperature 
systems, however, one should not use the true temperature the liquids 
and Kirkwood 1955; Oppenheim 1954). This important point since 
low-temperature liquids, such He* and increases decreases. For 

The final results are: 


Case 


Tes) 1/2 
1 


0 
Case 


where number molecules per cubic centimeter and the nuclear 
gyromagnetic ratio. have used the fact that the correlation times for 
self-diffusion are short compared with the Larmor period. The anisotropic 
interactions between molecules have been expanded 


¢(R) = da im(R) Ym (8, ¢); 


denotes Bessel function order the mean lifetime 

The best technique for measuring the pulse method Hahn (1950). 
This method has the additional advantage that enables direct and reliable 
measurement and hence direct comparison theory with experiment 
possible. Unfortunately, however, although has been measured for 
liquid (Bloom 1957) measurements were made. 
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Fig. have plotted the temperature dependence the diffusion 
constant liquid required explain the observed temperature dependence 
(Bloom 1957) according equation (4). this plot have neglected 
the temperature dependence introduced the integral; the integral varies 


| 


Fic. Plot predicted values diffusion constant for liquid function tem- 
perature using measured values and The diffusion constant plotted arbitrary 


units normalized the value 20° 


less than between 20° and 14° for liquid equilibrium with 
its vapor. The reason for this small effect that the temperature drops, 
and the density both increase and their effects g(R) tend cancel 
each other (Kirkwood, Lewinson, and Alder 1952). This term will, however, 
very important calculating the strong dependence ortho- 
concentration, the small difference the temperature dependence 
samples having different concentrations (Bloom 1957) and 
the effect pressure Results calculations bearing these experi- 
ments will presented subsequent publication. The reasonable tempera- 
ture dependence required explain the experimental behavior 
according this theory contrasted with predictions the 
Purcell-Pound (BPP) theory for liquids which gives where the 
viscosity. This expression predicts decreasing with decreasing temperature 
instead the observed increasing similar temperature dependence 
ours was obtained Moriya (1957) using different approach. 

not all surprising that the BPP theory breaks down completely 
low temperatures for case type relaxation. Their model sphere moving 
viscous medium used calculating correlation functions for molecular 
reorientation may fair approximation only for high values when 
the reorientational motion may described random walk process. For 
liquid such however, where only the state excited, the BPP 
model breaks down, single transition from one state another causes 
appreciable reorientation. Actually, all have done case extend 
the treatment gas the liquid case. 


~ 

| 

> 


NOTES 


The results for liquid equilibrium with its vapor are almost the 
same the BPP theory since, again, the integral equation (3) almost 
temperature independent. would therefore expect decrease with 
decreasing temperature, according measurements for (Garwin 
and Reich 1959). expected, however, that the statistics 
the atoms should become important and should cause increase 
strongly low temperatures upon the onset degeneracy. postpone 
discussion this difficult point since are here only concerned with 
dynamical effects. 

Finally, may remarked that experiments which would provide 
strong test this theory are the measurement function tempera- 
ture and pressure simple liquids such xenon and neon, where effects 
statistics are small. For these cases equation (3) should hold, and all the 
parameters appearing there are already known. 
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LETTERS THE EDITOR 


Under this heading brief reports important discoveries physics may published. These 
reports should not exceed 600 words and, for any issue, should submitted not later than six weeks 
previous the first day the month issue. proof will sent the authors. 


Frequency Measurement Standard Frequency Transmissions Against 
Caesium-Beam Resonator Standard* 


Since July 1958 caesium-beam frequency resonator has been operation the National 
Research Council laboratories (Kalra, Bailey, and Daams 1958, 1959). The local Essen ring 
quartz crystal oscillators (Kalra, Pattenson, and Thomson 1959) are calibrated periodically 
against this resonator standard. Since the day-to-day stability these oscillators the order 
these oscillators not need continuous monitoring against the caesium-beam 
frequency standard. equation the type computed each month for each 
oscillator, where the frequency any day, constant, and time days. 

Kalra, Pattenson, and Thomson (1959) have briefly noted the measurement the carrier 
frequency some the very low frequency stations with precision-controlled 
quencies phase comparison method. These measurements are made against one the 


TABLE I 


GBR 
16 ke/s 
August MSF, 6-hour 24-hour KK2XEI, 
1959 60 ke/s average* average 60 ke/s 
1 — 168 —167 N.M. 
2 — 167 — 173 N.M. 
4 — 162 N.M. — 86 
Hy — 158 — 169 — 83 
5 —161 — 163 — 87 
6 — 160 —157 — 87 
7 — 158 —161 —90 
— 161 —167 N.M. 
— 162 N.M. 
10 —155 — 157 
11 — 185 —170 —85 
12 N.M. —184 —85 
13 — 182 — 185 —179 —79 
14 N.M. — 176 — 83 
15 172 76 -177 N.M 
16 18) —176 -179 N.M 
17 177 ~ 169 
15 -173 —172 - 172 
1Y 175 175 171 ~ 81 
20 176 ~174 ~176 
21 177 17 -178 
22 172 178 ~174 N.M. 
23 176, 17% N.M. 
24 18) -178 82 
25 182 178 78 
26 180 175 N.M. 79 
27 182 174 N.M 7 
26 170 ~17% 174 80 
24 166 17% 169 N.M. 
20 172 170 170 N.M. 
31 172 166 170 76 
Mi ponthily 
mean —171 176 82 


Midmont bly 


mean ol WWV 104 


Nort: Measurements are made at Ottawa, Canada, using N.R.C, caesium 
beam frequency resonator as reference standard (with an assumed frequency ot 
9 192 631 770 cops). Frequency deviations from nominal are quoted in parts pet 
10! A negative sign indicates that the frequency is below nominal, 

N.M.: No measurement 


*Time of observations: 0100 to 0400 and 1000 to 1400 UT, 


*Jssued as N. B.C. No. 5498 


‘ 
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Essen ring oscillators. means the computed equation above, these measurements are 
later referred the caesium-beam frequency resonator standard. The time signals from some 
the standard frequency time stations are also monitored against these oscillators working 
clocks. The time signals and the carrier frequencies these stations are usually controlled 
the same crystal oscillator. The carrier frequency these stations thus computed with 
averaging period calendar month. For both techniques the over-all accuracy the 
frequency averaged around midmonth estimated Full details the techniques 
will published later. 

The accompanying table shows the frequency MSF (Rugby, England) GBR 
(Rugby, England) KK2XEI (Boulder, Colorado, U.S.A.) measured the 
phase technique; and the midmonthly mean frequency (Beltsville, Maryland, U.S.A.) 
and Mc/s, measured against the N.R.C. caesium-beam resonator standard. Values are 
quoted fractional parts 10", Af/fX where the departure from nominal and 
the nominal frequency. The reference frequency scale defined such that the N.R.C. caesium- 
beam resonator has frequency 192 631 770 c.p.s. negative sign indicates that frequency 
below nominal. 

Arrangements have been made for the monthly publication this Journal the received 
frequencies the above-mentioned signals measured this laboratory. 


Kara, S. N., Baitey, R., and Daams, H. 1958. Can. J. Phys. 36, 1442. 
1959. Nature (London), 183, 575. 
Kara, S. N., Pattenson, C. F., and THomson, M. M. 1959. Can. J. Phys. 37, 10. 
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